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Abstract 


Part  A  of  Phase  I  of  the  contract  was  to  develop  the  instrumentation  so  that  a 
dynamic  biomechanical  model  could  be  developed.  The  specific  objectives  were:  a)  to 
develop  Fastrak™  software  for  relative  pack-person  motion,  b)  construct  smaller  strap 
force  gauges,  c)  build  and  calibrate  a  moment  of  inertia  platform,  d)  modify  the  Load 
Carriage  compliance  tester  to  automate  two  degrees  of  freedom  and  e)  create  a  full  body 
mapping  of  mannequin.  This  report  describes  the  purposes  and  outputs  available  from  the 
pack-person  motion  (Section  2.0),  describes  the  development,  construction,  calibration 
and  protocol  for  use  of  smaller  strap  sensors  (Section  3.0)  and  the  moment  of  inertia 
platform  (Section  4.0),  and  the  development  and  steps  involved  in  modifying  the  LC 
compliance  tester  (Section  5.0)  and  the  mannequin  mapping  (Section  6.0).  For  the  most 
part,  tasks  were  developmental  and  construction-based  with  no  data  analyses,  other  than 
to  confirm  the  accuracy  and  precision  of  the  instrumentation. 

This  report  deals  only  with  Part  A  consisting  of  five  sub-parts  of  the  contract.  Part 
B  involved  changing  the  technical  manuals  based  on  the  upgrades  stated  within  this 
report.  Part  C  of  this  contract  is  under  separate  cover  and  was  to  develop  a  long  range 
plan  and  budget  for  dynamic  biomechanical  modeling.  Part  D  was  to  assist  with  the 
NATO  HFM  Specialist  Meeting  entitled  “Soldier  Mobility:  Innovations  in  Load  Carriage 
System  Design  and  Evaluation”  held  on  27-29  June  2000.  Part  D  is  described  by  the 
NATO  RTO  MP-56  Technical  Proceedings  Report  entitled:  “Soldier  Mobility: 
Innovations  in  Load  Carriage  System  Design  and  Evaluation”. 
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Resume 


La  partie  A  de  la  phase  I  du  contrat  avait  pour  but  la  mise  au  point  d’instruments 
a  fin  de  permettre  l’elaboration  d’un  modele  biomecanique  dynamique.  Elle  visait  les 
objectifs  precis  suivants  :  a)  developper  le  logiciel  Fastrakmc  alin  de  determiner  le 
mouvement  relatif  entre  un  sac  a  dos  et  une  personne,  b)  construire  de  plus  petits 
dynamometres  de  sangle,  c)  fabriquer  et  etalonner  une  plate-forme  de  moment  d’inertie, 
d)  modifier  l’appareil  d’essai  de  conformite  de  transport  tactique  afin  d’automatiser 
2  degres  de  liberte  et  e)  creer  une  representation  integrale  du  corps  du  mannequin.  Le 
present  rapport  decrit  la  raison  d’etre  du  logiciel  et  de  ses  sorties  relativement  au 
mouvement  entre  le  sac  a  dos  et  la  personne  (section  2.0);  la  mise  au  point,  la  fabrication, 
l’etalonnage  et  le  protocole  d’utilisation  de  plus  petits  capteurs  de  sangle  (section  3.0) 
ainsi  que  de  la  plate-forme  de  moment  d’inertie  (section  4.0);  le  developpement  de 
l’appareil  d’essai  de  conformite  de  transport  tactique  et  les  etapes  suivies  pour  le  modifier 
(section  5.0)  et  la  representation  du  mannequin  (section  6.0).  La  majeure  partie  du  travail 
consistait  en  des  taches  de  developpement  et  de  fabrication  et  ne  necessitait  aucune 
analyse  de  donnees,  autre  que  pour  etablir  la  justesse  et  la  precision  des  instruments. 

Le  present  rapport  traite  seulement  de  la  partie  A  du  contrat,  qui  comprend  cinq 
sous-parties.  La  partie  B  consiste  a  modifier  les  manuels  techniques  en  tenant  compte  des 
mises  a  niveau  indiquees  dans  le  present  rapport.  La  partie  C  du  contrat  est  fournie 
separement  et  porte  sur  l’etablissement  d’un  plan  et  d’un  budget  a  long  terme  pour 
l’elaboration  d’un  modele  biomecanique  dynamique.  La  partie  D  a  ete  redigee  pour  la 
reunion  des  specialistes  HFM  de  l’OTAN,  qui  est  intitulee  Soldier  Mobility:  Innovations 
in  Load  Carriage  System  Design  and  Evaluation,  et  qui  a  eu  lieu  du  27  au  29  juin  2000. 
Une  description  en  est  donnee  dans  le  compte  rendu  technique  ORT  MP-56  intitule  «  La 
mobilite  du  combattant :  innovations  dans  la  conception  et  l’evaluation  des  gilets 
d’ intervention  [document  en  anglais  seulement]  ». 
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Executive  Summary 


The  purpose  of  this  contract  was  to  upgrade  the  equipment  and  improve  the 
instrumentation  in  support  of  research  and  development  into  dynamic  biomechanical 
modeling  for  better  LC  design  and  system  evaluation,  and  development  of  a  new 
STANAG  for  load  carriage  and  is  reported  in  two  phases.  The  Phase  1  report  deals 
primarily  with  hardware  components  needed  to  develop  a  dynamic  biomechanical  model. 
Each  subsection  of  the  study  will  be  described  in  summary  below. 

A.l  Fastrak™  Software  for  Relative  Pack-Person  Motion 

The  Fastrak™  is  a  six-degree  of  freedom  tracking  instrument  that  uses 
electromagnetic  fields  to  determine  the  position  and  orientation  of  a  remote  object.  Our 
original  set-up  had  the  source  mounted  on  the  mannequin  and  either:  1)  one  sensor  placed 
in  the  pack  for  pack-person  motion  or  2)  four  sensors  placed  in  the  soldier’s  kit  for  kit- 
person  motion.  This  set-up  provided  3D  pack-person  motion  only;  however,  this  set-up 
did  not  tell  us  whether  the  pack  and  kit  were  in-synch  with  the  person.  Hence  it  was 
necessary  to  revise  the  dynamic  biomechanical  model.  The  set-up  and  software 
development  was  revised  with  the  addition  of  a  second  sensor  as  a  ground  reference  point. 
This  approach  allowed  us  to  detennine  if  the  pack  and  person  were  in-phase  or  out-of¬ 
phase,  and  the  relative  amplitudes  of  their  motions.  A  physical  structure  was  mounted  in 
the  ceiling  to  hold  the  sensor.  Most  of  the  work,  however,  dealt  with  the  revision  of  the 
software.  The  new  software  analysis  program  can  be  used  for  several  purposes:  a)  it  can 
examine  the  frequency  content  of  the  pack  and  person  during  standardized  motions  on  the 
LC  Simulator,  b)  it  can  be  used  to  detennine  whether  a  suspension  system  can  be 
designed  to  transfer  mechanical  energy  within  the  gait  cycle  to  assist  the  person  in 
walking,  and  c)  it  can  be  used  to  create  a  better  dynamic  biomechanical  model  (DBM) 
because  it  can  help  model  the  actions  of  the  pack-person  interface  force  requirements. 

A.2  Construct  Small  Strap  Force  Gauges 

The  current  design  for  strap  force  sensors  are  dog-bone  shaped  and  made  out  of  an 
aluminium  base  plate  with  rosette  strain  gauges  that  require  a  minimum  of  6  cm  of  strap 
length.  These  strap  force  sensors  are  linear  and  robust  and  the  design  works  well  for 
shoulder  straps  and  the  waist  belt.  However,  they  are  not  suitable  for  load-lifter  straps 
and  hip  stabilizer  straps  because  of  length  of  strap  available  and  the  radius  of  curvature  at 
some  locations  around  the  shoulders  and  hips.  The  purpose  of  this  phase  was  to  construct 
and  calibrate  smaller  strap  tension  sensors  to  access  the  difficult-to-measure  straps  (i.e., 
load  lifters,  hip  stabilizers).  The  goal  was  to  develop  eight  strap  gauges  that  could 
withstand  the  same  tension  ranges  (~  5 ON  to  100N)  but  be  used  in  tight  locations.  These 
sensors  were  constructed,  calibrated  and  checked  for  linearity  and  DC  offset.  Physical 
length  of  the  new  gauges  is  18  mm  and  they  can  be  sewn  into  position.  Each  sensor  is 
now  available  for  use. 
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A.3  Construction  of  a  Moment  of  Inertia  Platform 


One  requirement  of  dynamic  biomechanical  modeling  (DBM)  is  to  know  the 
location  of  the  centre  of  gravity  and  the  moment  of  inertia  of  a  pack.  Due  to  the  non- 
uniform  nature  of  a  backpack,  pinpointing  the  centre  of  gravity  can  be  difficult  and 
mathematically  calculating  the  moment  of  inertia  of  a  non-uniformly  shaped  and  weighted 
object  requires  complex  equations.  A  trifilar  pendulum  provides  a  relatively  simple  and 
quick  method  of  locating  the  centre  of  gravity  and  the  moment  of  inertia  of  non- 
homogenous  objects.  A  trifilar  pendulum  was  built  which  consisted  of  a  round  laminated 
plywood  circular  plate  that  was  suspended  in  the  horizontal  plane  by  three  cables  of  equal 
length  equidistantly  placed  from  the  centre  of  the  plate.  When  the  pendulum  is  twisted 
around  a  vertical  axis  running  through  the  centre  of  the  plate,  it  oscillates  with  a  period 
dependent  on  the  following  factors:  the  length  of  the  cables,  the  radius  of  the  plate,  the 
mass  of  the  plate,  the  local  force  of  gravity,  and  the  moment  of  inertia  of  the  plate.  This 
trifilar  pendulum  was  calibrated  and  will  be  used  to  find  the  centre  of  gravity  and  moment 
of  inertia  of  backpacks  for  the  dynamic  biomechanical  model. 

A.4  Modify  Load  Carriage  Compliance  Tester 

The  compliance  tester  was  originally  developed  to  assess  the  stiffness  of  equipment 
worn  on  the  torso  using  quasi-static  motion.  Movement  of  the  test  torso  was  achieved 
using  a  series  of  hand  operated  pulleys  and  cables.  This  required  an  operator  to  attempt 
to  produce  a  constant  rate  of  displacement  that  introduced  some  variability  in  the 
methodology  and  confined  testing  to  quasi-static  conditions.  The  purpose  for  the  revised 
LC  Compliance  tester  design  modification  was  to  support  the  development  of  a 
biomechanical  model  of  human  load  carriage.  The  device  will  furnish  input  values  for  the 
dynamic  stiffness  variables  during  initial  modelling  stages.  At  later  stages  of  the  model 
development,  the  apparatus  can  be  used  to  validate  the  model  by  comparing  the  model 
predicted  stiffness  responses  to  a  different  range  of  motion  profiles.  Only  2D  changes 
were  made  to  the  LC  Compliance  tester  at  this  time  with  the  3D  function  to  be  completed 
later.  These  changes  to  the  automated  programmable  motion  control  system  create  a 
number  of  functions  not  possible  with  the  previous  system.  These  functions  include:  1) 
highly  repeatable  motion  profiles  independent  of  the  operator;  2)  determination  of  the 
system  stiffness  under  dynamic  conditions;  3)  provision  for  quantitative  validation  of 
models  of  load  carriage  devices;  4)  determination  of  the  frequency  response  of  LC 
suspension  systems;  and  5)  creation  of  an  automated  test  cell  requiring  minimal  operator 
expertise  and  low  cost  for  possible  sale  to  support  other  countries’  modeling  efforts. 

A.5  Create  Full  Body  Mapping  of  Mannequin 

Our  previous  work  with  static  biomechanical  modelling  yielded  relatively  good 
results  if  the  pack  was  kept  simple  (i.e.  shoulder  straps  and  initial  waist  belt  designs). 
However,  as  more  pack  elements  were  added  (i.e.  load  lifter  straps,  hip  stabilizers  etc)  the 
models  were  less  effective,  primarily  because  of  the  difficulty  in  properly  modeling  the 
pack-person  interface  forces.  The  purpose  of  this  equipment  upgrade  was  to  allow  a 
pressure  measurement  system  to  be  used  as  an  input  to  the  dynamic  biomechanical  model. 
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This  was  accomplished  in  two  ways:  1)  write  new  software  to  allow  direct  access  to  the 
raw  data  from  the  TeKScan™  thus  creating  the  potential  to  extract  specific  contact  force 
histories  from  Tekscan™  and  2)  to  develop  the  ability  to  detennine  nonnal  force  vectors 
and  locate  them  in  space  on  the  body  by  conducting  a  total  body  mapping  of  the  50th 
percentile  male  mannequin.  For  the  first  objective,  a  combination  of  MatLab  (a  numerical 
computation  and  data  visualization  program)  and  Microsoft  Excel  macros  was  written. 
This  software-based  approach  was  developed  to  manipulate  the  files  into  the  appropriate 
formats  compatible  with  MSC  Visual  Nastran™  4D  and  2D  Working  Model.  A  full  body 
mapping  of  the  50th  percentile  mannequin  is  necessary  to  locate  the  senseis  in  3D  space, 
and  it  will  allow  the  calculation  of  the  normal  reaction  force.  The  LC  mannequin  was 
scanned  so  that  the  full  body  mapping  files  are  now  accessible  by  Visual  Nastran  4D  for 
validation  of  the  waist  belt  model.  The  resolution  of  this  scan  was  2mm. 
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Sommaire 


Le  contrat  visait  a  mettre  a  niveau  l’equipement  et  a  ameliorer  les  instruments,  qui 
serviront  a  la  R  et  D  sur  la  modelisation  biomecanique  dynamique  en  vue  d’ameliorer  la 
conception  et  1’ evaluation  des  systemes  de  transport  de  charge,  et  a  e  laborer  un  nouveau 
STANAG  sur  le  transport  de  charge.  II  comporte  deux  phases.  Le  rapport  sur  la  phase  1 
porte  principalement  sur  le  materiel  necessaire  a  T  elaboration  d’un  modele  biomecanique 
dynamique.  Une  breve  description  de  chaque  sous-section  de  T  etude  est  donnee  ci- 
dessous. 

A.l  Logiciel  FastrakMC  pour  mesurer  le  mouvement  relatif  entre  un  sac  a  dos  et 
une  personne 

FastrakMC  est  un  instrument  de  localisation  a  six  degres  de  liberte  qui  utilise  des 
champs  electromagnetiques  pour  determiner  la  position  et  l’orientation  d’un  objet  eloigne. 
Dans  la  configuration  originale,  la  source  etait  installee  sur  le  mannequin  et  1)  un  capteur 
etait  place  dans  le  sac  a  dos  pour  detecter  le  mouvement  sac/personne  ou  2)  quatre 
capteurs  etaient  places  dans  le  fourbi  du  soldat  pour  detecter  le  mouvement 
fourbi/personne.  Cette  configuration  permettait  de  connaitre  seulement  le  mouvement  sac 
a  dos/personne  en  trois  dimensions;  toutefois,  elle  ne  permettait  pas  de  savoir  si  le 
mouvement  du  sac  et  du  fourbi  etait  synchronise  avec  celui  de  la  personne.  Done,  il  a 
fallu  modifier  le  modele  biomecanique  dynamique.  La  configuration  et  le  logiciel  ont  ete 
modifies  en  tenant  compte  de  l’ajout  d’un  deuxieme  capteur  comme  point  de  reference  au 
sol.  Cette  methode  a  permis  de  determiner  si  les  mouvements  du  sac  a  dos  et  de  la 
personne  etaient  en  phase  ou  dephases  et  de  connaitre  1’ amplitude  relative  de  ces 
mouvements.  Une  structure  a  ete  posee  au  plafond  pour  contenir  le  capteur.  La  majeure 
partie  du  travail  a  ete  consacree  a  la  modification  du  logiciel.  Le  nouveau  logiciel 
d’ analyse  peut  etre  utilise  a  plusieurs  fins,  soit  a)  examiner  les  frequences  produites  par  le 
sac  et  la  personne  durant  des  mouvements  normalises  effectues  avec  le  simulateur  LC, 

b)  determiner  si  un  systeme  de  suspension  capable  de  transferer  l’energie  mecanique  au 
cycle  de  marche  peut  etre  concu  afin  de  faciliter  les  deplacements  de  la  personne  et 

c)  ameliorer  le  modele  biomecanique  dynamique  (DBM)  en  modelisant  les  forces 
exercees  au  niveau  de  l’interface  sac-personne. 

A.2  Fabrication  de  dynamometres  de  sangle  plus  petits 

Actuellement,  les  capteurs  de  force  de  sangle  sont  en  forme  d’os  a  chien  et 
consistent  en  une  plaque  d’aluminium  dotee  d’extensometres  en  rosette  qui  necessitent 
une  sangle  d’au  mo  ins  6  cm  de  longueur.  Ces  capteurs  de  force  de  sangle  sont  lineaires  et 
robustes,  et  leur  conception  est  bien  adaptee  aux  sangles  portees  aux  epaules  et  a  la  taille. 
Toutefois,  ils  ne  conviennent  pas  aux  sangles  de  levage  de  charge  et  aux  sangles  de 
stabilisation  aux  hanches  en  raison  de  la  longueur  de  sangle  disponible  et  du  rayon  de 
courbure  a  certaines  parties  des  epaules  et  des  hanches.  Cette  phase  visait  a  fabriquer  et  a 
etalonner  des  capteurs  de  traction  de  sangle  plus  petits  afin  de  permettre  la  mesure  sur  les 
sangles  difficiles  d’acces  (c.-a-d.  les  sangles  de  levage  de  charge  et  de  stabilisation  aux 
hanches).  L’objectif  etait  de  mettre  au  point  huit  dynamometres  de  sangle  qui  peuvent 
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tenir  a  la  meme  plage  de  traction  (~  50  N  a  100  N)  tout  en  etant  utilises  dans  des  endroits 
exigus.  Une  fois  les  capteurs  fabriques  et  etalonnes,  on  a  verifie  la  linearite  et  le  decalage 
c.c.  Les  nouveaux  capteurs  mesurent  18  mm  de  long  et  peuvent  etre  cousus  en  place.  Tous 
les  capteurs  sont  maintenant  prets  a  etre  utilises. 

A.3  Fabrication  d’une  plate-forme  de  moment  d’inertie 

Le  modele  biomecanique  dynamique  (DBM)  exige  que  l’on  connaisse  le  centre  de 
gravite  et  le  moment  d’inertie  d’un  sac  a  dos.  En  raison  du  caractere  asymetrique  d’un  sac 
a  dos,  il  est  difficile  de  determiner  avec  precision  le  centre  de  gravite  de  ce  dernier  et, 
pour  calculer  le  moment  d’inertie  d’un  objet  dont  la  forme  et  le  poids  sont  inegaux,  il  faut 
recourir  a  des  equations  complexes.  Un  pendule  trililaire  constitue  un  moyen  relativement 
simple  et  rapide  pour  determiner  le  centre  de  gravite  et  le  moment  d’inertie  d’objets  non 
homogenes.  On  en  a  fabrique  un  en  utilisant  une  plaque  circulaire  en  lamelles  de 
contreplaque  et  on  l’a  suspendu  sur  le  plan  horizontal  en  utilisant  trois  cables  de  meme 
longueur  places  a  distance  egale  du  centre  de  la  plaque.  Lorsque  le  pendule  tourne  autour 
d’un  axe  vertical  passant  par  le  centre  de  la  plaque,  il  oscille  selon  une  periode  fondee  sur 
les  facteurs  suivants  :  la  longueur  des  cables,  le  rayon  de  la  plaque,  la  masse  de  la  plaque, 
la  force  de  gravite  locale  et  le  moment  d’inertie  de  la  plaque.  Le  pendule  trililaire  a  ete 
etalonne  et  sera  utilise  pour  determiner  le  centre  de  gravite  et  le  moment  d’inertie  des  sacs 
a  dos  aux  fins  du  modele  biomecanique  dynamique. 

A.4  Modifier  l’appareil  d’essai  de  conformite  de  transport  tactique 

A  l’origine,  l’appareil  d’essai  de  conformite  a  ete  developpe  pour  evaluer  la  rigidite  de 
l’equipement  porte  sur  le  torse  dans  des  conditions  de  mouvement  quasi-statiques.  Le 
torse  d’essai  etait  mu  au  moyen  d’une  serie  de  poulies  et  de  cables  manipules  par 
l’operateur.  Ce  dernier  devait  alors  essayer  de  produire  une  vitesse  de  deplacement 
constante.  Cela  a  eu  pour  effet  de  creer  une  certaine  variation  dans  la  methode  et  de 
limiter  les  essais  a  des  conditions  quasi-statiques.  Les  modifications  apportees  a  l’appareil 
d’essai  de  conformite  de  transport  tactique  avaient  pour  but  d’appuyer  l’elaboration  d’un 
modele  biomecanique  de  transport  de  charge  par  une  personne.  Durant  les  premieres 
etapes  de  la  modelisation,  l’appareil  fournira  les  valeurs  d’entree  pour  les  variables  de 
rigidite  dynamique.  Durant  des  etapes  ulterieures,  il  pourra  etre  utilise  pour  valider  le 
modele  en  comparant  les  reponses  de  rigidite  prevues  avec  differents  profils  d’ amplitude 
de  mouvement.  A  l’heure  actuelle,  seules  les  modifications  bidimensionnelles  ont  ete 
apportees  a  l’appareil;  la  fonction  tridimensionnelle  sera  ajoutee  plus  tard.  Les 
changements  apportes  au  systeme  automatise  programmable  de  commande  des 
mouvements  ont  donne  acces  a  des  fonctions  qui  n’existaient  pas  dans  le  systeme 
precedent,  notamment :  1)  profils  de  mouvements  tres  reproductibles,  independants  de 
l’operateur,  2)  determination  de  la  rigidite  du  systeme  dans  des  conditions  dynamiques, 

3)  validation  quantitative  de  modeles  de  systemes  de  transport  de  charge,  4)  determination 
de  la  reponse  en  frequence  des  systemes  de  suspension  de  transport  de  charge  et 
5)  creation  d’une  cellule  d’essai  automatisee  peu  couteuse  ne  necessitant  pas  beaucoup  de 
competences  de  la  part  de  l’operateur,  qui  peut  etre  vendue  afin  de  soutenir  les  activites 
de  modelisation  d’autres  pays. 
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A.5  Creer  une  representation  integrate  du  corps  du  mannequin 

Nos  recherches  anterieures  relatives  au  modele  biomecanique  statique  ont  donne 
des  resultats  relativement  bons  lorsque  le  sac  a  dos  etait  simple  (c.-a-d.  modeles  initiaux 
avec  sangles  d’epaule  et  ceinture).  Cependant,  avec  l’ajout  d’elements  supplementaires  au 
sac  (sangles  de  levage  de  charge,  sangles  de  stabilisation  aux  hanches,  etc.),  les  modeles 
ont  perdu  de  leur  efficacite,  principalement  en  raison  de  la  difficulty  a  modeliser 
correctement  les  forces  a  l’interface  sac-personne.  La  mise  a  niveau  de  l’equipement  avait 
pour  but  de  pennettre  l’utilisation  d’un  systeme  de  mesure  de  la  pression  cornme  entree 
du  modele  biomecanique  dynamique.  A  cette  fin,  les  deux  modifications  suivantes  ont  ete 
apportees  :  1)  creation  d’un  nouveau  programme  pour  pennettre  faeces  direct  aux 
donnees  brutes  de  l’appareil  TekscanMC,  ce  qui  permet  d’extraire  des  historiques 
specifiques  des  forces  de  contact  a  partir  du  TekscanMC  et  2)  creation  d’une  representation 
integrale  du  corps  d’un  mannequin  representant  le  50e  centile  de  sexe  masculin  afin  de 
determiner  les  vecteurs  de  force  normaux  et  de  les  localiser  sur  le  corps.  Pour  atteindre  le 
premier  objectif,  des  macros  ont  ete  ecrites  en  MatLab  (programme  de  calcul  numerique 
et  de  visualisation  de  donnees)  et  en  Excel,  de  Microsoft.  Cette  methode  fondee  sur 
logiciel  a  ete  elaboree  afin  de  convertir  les  fichiers  aux  formats  appropries,  compatibles 
avec  le  modele  de  travail  2D  et  4D  MSC  Visual  NastranMC.  Une  representation  integrale 
du  corps  du  mannequin  male  du  50e  centile  est  necessaire  pour  localiser  les  capteurs  a 
cellules  dans  un  espace  tridimensionnel  afin  de  pennettre  le  calcul  de  la  force  de  reaction 
normale.  Le  balayage  du  mannequin  utilise  pour  le  transport  des  charges  a  ete  effectue,  et 
les  fichiers  de  mappage  du  corps  entier  sont  maintenant  accessibles  au  moyen  de  Visual 
Nastran  4D  pour  fins  de  validation  du  modele  de  ceinture.  La  resolution  du  balayage  est 
de  2  mm. 
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1.0  Introduction 


The  work  on  the  Dynamic  Biomechanical  Model  (DBM)  has  been  underway  after 
the  duties  related  to  the  NATO  Load  Carriage  Workshop  were  completed.  One  of  the 
outcomes  from  the  NATO  HFM  Specialist  Meeting  entitled  “Soldier  Mobility: 
Innovations  in  Load  Carriage  System  Design  and  Evaluation”  held  on  27-29  June  2000, 
was  to  work  toward  the  development  of  a  NATO  STANAG  based  on  the  new  knowledge 
developed  by  various  countries  involved  in  research  related  to  load  carriage.  The 
contribution  Canada  can  make  toward  this  aim  is  the  development  of  the  minimal 
biomechanical  human  tolerance  limits  based  on  such  variables  as  shoulder  reaction  forces, 
lumbar  reaction  forces  and  pressure  tolerance  limits.  One  method  to  acquire  these  data  is 
to  develop  a  dynamic  biomechanical  model  that  can  provide  direct  input  of  load  control 
and  load  transfer  conditions  from  soldier  trials  with  different  systems.  Eventually  this 
biomechanical  model  could  be  used  in  conjunction  with  a  portable  measurement  system 
to  provide  on-line  feedback  about  these  variables  to  identify  the  magnitudes  at  which 
local  muscular  fatigue  or  pain  occurs  that  can  affect  operational  effectiveness.  Using  this 
modified  equipment  and  a  DBM  design  tool,  there  is  also  the  potential  to  enhance  our 
understanding  and  thus  optimization  of  the  transfer  of  mechanical  energy  between  the 
pack  and  person  in  the  design  of  new  systems  thus  creating  “a  smart  LC  system”. 


1.1  Purpose 

The  purpose  of  this  part  of  the  study  was  to  upgrade  the  equipment  and 
instrumentation  in  support  of  research  and  development  into  dynamic  biomechanical 
modeling  for  better  design  and  LC  system  evaluation,  and  development  of  a  new 
STANAG  for  load  carriage.  Only  the  equipment  prepared  in  Phase  1  Part  A  of  this  multi¬ 
phase  study  will  be  described  in  this  report,  namely: 

1.0  Introduction 

1.0  Fastrak™  software  for  relative  pack-person  motion 
2.0  Construct  small  strap  force  gauges 
3.0  Construction  of  a  moment  of  inertia  platform 
4.0  Modify  Load  Carriage  compliance  tester 
5.0  Create  full  body  mapping  of  mannequin. 
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2.0  Fastrak™  software  for  relative  pack-person  motion 

2.1  Background 


In  the  previous  research  and  development  phases  to  develop  the  LC  Simulator 
(Stevenson  et  ah,  1995;  1997),  a  four  channel  electromagnetic  sensing  system  was  used  to 
detect  relative  pack-person  motion.  The  FASTRAK™,  by  Polhemus  Navigation 
Incorporated,  Colchester,  Vennont,  USA  is  a  six-degree  of  freedom  tracking  instrument 
that  uses  electromagnetic  fields  to  detennine  the  position  and  orientation  of  a  remote 
object.  The  system  consists  of  three  components:  a  systems  electronic  unit  (SEU),  a 
source,  and  up  to  four  sensors.  The  source  generates  a  low  frequency  electromagnetic 
field  that  is  detected  by  the  sensors.  The  SEU  calculates  the  three  dimensional  position 
and  orientation  of  the  sensor  in  this  field  relative  to  the  source.  This  is  accomplished  by 
interpreting  the  interaction  of  the  electromagnetic  fields  between  three  sets  of  orthogonal 
coils  contained  in  both  the  source  and  the  sensor  (An  et  ah,  1988;  Day  et  ah,  1998; 
Polhemus,  1995).  The  sensors  can  be  detected  by  the  source  up  to  a  maximum  separation 
distance  of  305  cm,  but  it  is  recommended  that  all  testing  occur  within  76  cm  of  the 
source.  The  specifications  state  that  any  metal  within  the  76  cm  envelope  may  cause 
inaccurate  position  and  orientation  data,  but  calibration  performed  by  direct  comparison 
of  Fastrak™  positional  data  with  data  collected  from  an  opto-electric  positional  recording 
system  (Optotrak™  by  Northern  Digital  Incorporated)  with  high  precision  (RMS  error 
<0.01  mm)  provided  an  RMS  error  for  Fastrak™  data  of  0.65  mm  (Stevenson  et  ah, 

1995).  Figure  2-1  is  a  picture  of  the  source  and  a  sensor. 


Figure  2-1.  Fastrak™  sensor  (left)  and  source  (centre).  Data  are  taken  via  the  RS232  port 
connection  (right)  to  a  portable  computer. 
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In  the  APLCS  (Advanced  Personal  Load  Carriage  System)  project,  the  FASTRAK™ 
source  was  affixed  with  nylon  screws  to  the  underside  of  the  left  arm  of  the  mannequin  so 
that  the  Z  axes  was  superimposed  on  the  mannequin’s  longitudinal  axis.  This  allowed  the 
source  to  represent  the  orientation  of  the  body  axis  as  the  global  reference  for  the  pack 
motion.  The  source  was  oriented  such  that  the  X  axis  corresponded  to  the  sagittal  plane, 
the  Y  axis  corresponded  to  the  transverse  plane,  and  the  Z  axis  corresponded  to  the  frontal 
plane.  In  other  words,  the  front/back  (anterior/posterior)  motions  were  recorded  by  the  ± 
X  direction,  side-to-side  (medial/lateral)  motions  by  the  ±  Y  direction  and  up-down 
(superior/inferior)  motions  were  represented  by  ±  Z  direction.  Pack  resultant  motion  was 
calculated  for  pack-person  relative  displacement. 


Figure  2-2  Original  set-up  for  assessment  of  relative  motion  for  Fastrak™  system. 

For  typical  pack-based  protocol  on  the  LC  Simulator,  the  source  is  mounted 
securely  to  the  left  ann  with  angled  plexiglass  and  plastic  screws.  A  Fastrak™  sensor 
would  be  placed  securely  into  the  superior  polystyrene  surface  of  the  LC  System  payload 
in  order  to  acquire  displacement  data  for  the  payload  with  respect  to  the  source.  Data 
would  be  collected  for  10  seconds  at  55  Hz  on  300  second  intervals  over  a  15  minute  test 
duration.  Coordinate  transformation  was  applied  to  translate  from  the  superior  payload 
sensor  location  to  the  loaded  pack  centre  of  gravity  to  estimate  resultant  displacement 
vectors  between  the  centre  of  gravity  of  the  loaded  LC  system  and  the  mannequin.  The 
typical  output  from  the  previous  data  analysis  is  shown  in  Figure  2-3. 
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Figure  2-3.  A  superior  LC  system  (top  10th  percentile)  would  have  relative  motion 
values  between  the  pack  and  person  of  1 .32  mm  in  the  x  direction,  1 .06  mm  in  the  y 
direction  and  8.16  mm  in  the  z  direction. 

Based  on  the  database  of  load  carriage  systems  tested,  the  amount  of  x,  z,  and  r 
(total  relative  motion)  was  related  to  soldier  hip  discomfort  scores  (Bryant  et  ah,  1997). 
Increased  displacement  of  the  LC  system,  with  respect  to  the  soldier,  can  also  decrease 
agility  due  to  reduced  load  control.  Unrestrained  displacement  of  marching  order  kit  and 
LC  equipment  also  leads  to  stability  problems  and  local  discomfort  due  to  repeated 
collisions  between  the  soldier  and  the  pack  and  other  kit  items. 
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2.2  Rationale  for  Revised  Relative  Motion  Data 


Martin  et  al.  (1982)  felt  that  if  the  load  moved  independent  of  the  soldier  during 
load  carriage,  more  energy  must  be  used  to  control  the  pack  and  greater  fatigue  would 
develop.  Therefore,  they  hypothesized  that  control  of  the  relative  motion  between  pack 
and  person  was  important  so  that  load  carriage  system  moved  with,  rather  than  oppose  the 
wearer’s  nonnal  gait  pattern.  It  is  this  concept  that  lead  Queen’s  to  propose  relative 
pack-person  motion  that  represent  a  superior  top  10%  pack  (Table  2-1)  as  the  standard  of 
acceptability  for  future  load  carriage  systems.  In  other  word,  pack-person  motion  that 
was  as  small  as  possible  was  deemed  the  best  situation  for  load  control  and  load  transfer. 

Table  2-1.  Pack-Person  delta  displacement  values  that  represent  the  10th  ,  50th  and 
90th  percentile  of  seventeen  packs  from  database. 


Direction 

10th  (low) 
Decile 

Mean  (50th) 
Decile 

90th  (high) 
Decile 

X  (mm) 

1.32 

6.82 

12.33 

Y  (mm) 

1.06 

3.83 

66.60 

Z  (mm) 

7.47 

11.32 

15.17 

R  (mm) 

8.16 

14.06 

19.97 

Kram  in  1991  attempted  to  explain  why  Asian  people  carried  loads  using  compliant 
springy  bamboo  or  wooden  poles  on  the  shoulders  (Figure  2-3).  At  the  recent  NATO 
conference,  Kram  further  hypothesized  that  these  springy  carrying  poles  acted  like  a  spring- 
mass  system  where  the  poles  should  vibrate  at  1  Hz  if  the  person  is  walking  at  a  3  Hz  pace. 
This  would  create  a  natural  frequency  of  vibration  equal  to  the  square  root  of  k/m  where  k 
is  the  spring  stiffness  and  m  is  the  mass  supported  by  the  spring.  Although  Kram  (1991) 
was  not  able  to  demonstrate  his  theory  of  energetics  for  this  system,  he  felt  that  compliant 
suspension  systems  would  allow  better  energy  transfer  to  take  place. 


Figure  2-4  The  walker  causes  the  poles  to  vibrate  out  of  phase  this  each  step  so  there  is 
transfer  of  energy  from  the  pole  to  the  person  (with  permission  of  R.  Kram,  2000). 
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Kram  (2000)  described  the  relationship  between  gravitational  potential  energy 
(mgh)  and  kinetic  energy  (1/2  mv  )  as  an  inverted  pendulum-like  mechanism  during 
walking  whereby  the  interchange  saved  the  body  from  work  to  lift  and  accelerate  the 
center  of  mass  during  each  step  (Wcom).  Based  on  measuring  the  metabolic  cost  of  normal 
walking,  most  humans  appear  to  recover  about  65  %  of  the  energy  from  step  to  step,  less 
than  the  75%  Heglund  et  al.  (1995)  found  with  trained  head  load  carriers.  Heglund  et  al. 
(1995)  found  that  experienced  head  load  carriers  could  carry  up  to  20%  of  their  body 
mass  on  their  heads  with  no  measurable  increase  in  metabolic  energy  consumption  but 
Kram  (2000)  found  a  linear  increase  of  the  metabolic  relationship  for  loads  carried  at 
waist  level  over  the  same  20%  load  increase.  Both  researchers  found  that  recovery  of 
mechanical  energy  increased  with  heavier  loads  but  the  mechanical  work  performed  on 
the  center  of  mass  plus  load  did  not  increase  for  loads  greater  than  20%  of  body  mass. 
Kram  (2000)  hypothesized  that  the  metabolic  cost  of  walking  was  detennined  by  three 
factors:  cost  of  performing  external  work,  cost  of  generating  muscular  force  to  support 
body  weight  and  cost  of  swinging  the  limbs.  These  costs  can  change  with  tasks 
perfonned.  For  example  during  running,  tendons  can  return  as  much  as  93%  of  the 
energy  stored  in  them  (Alexander,  1988)  thus  saving  the  muscles  from  the  task  repeatedly 
perfonning  work  to  lift  and  accelerate  the  center  of  mass  during  each  step  thus 
dramatically  reduces  metabolic  cost. 

Kram  (2000)  used  these  concepts  to  describe  the  most  appropriate  suspension 
system  of  load  carrying  devices  similar  to  using  springy  bamboo  or  wooden  poles.  Using 
loads  on  the  shoulders,  Kram  (1991)  demonstrated  that  peak  vertical  ground  reaction 
force  were  around  2.5  times  body  weight  but  properly  tuned  carrying  pole  applied  a 
steady  and  moderate  force  with  force  fluctuations  less  than  10%,  however,  there  was  no 
metabolic  energy  savings  with  springy  poles.  Kram  challenged  designers  to  incorporate 
compliant  suspension  systems  into  load  carrying  systems  used  by  soldiers  and  measure 
the  impact  of  this  change. 

The  concept  of  a  compliant  suspension  system  was  part  of  the  impetus  to  develop 
the  lateral  rods  in  the  Clothe  the  Soldier  (CTS)  system  (Reid  et  al.,  2000).  Based  on  the 
expert  opinion  of  designer  Bill  Ostrum  from  Ostrum  Outdoors,  he  recognized  that  the 
rods  “worked  with”  the  trekker  and  seemed  to  even  out  the  energy  requirements.  We 
could  not  evaluate  this  potential  at  Queen’s  because  we  did  not  have  the  appropriate 
protocol  or  software  to  determine  the  natural  frequency  of  the  pack  and  person.  However, 
Queen’s  researchers  were  able  to  show  that  rods  under  tension  could  change  the  required 
extensor  force  by  the  erector  spinae  and  shift  some  of  the  force  requirements  from  the 
shoulders  to  the  waist. 


2.3  Purpose  of  the  Fastrak™  Software 

This  is  a  new  software  analysis  program  which  can  be  used  for  several  purposes: 
a)  it  can  examine  the  frequency  content  of  the  pack  and  person  during  standardized 
motions  on  the  LC  Simulator,  b)  it  can  be  used  to  determine  whether  a  suspension  system 
can  be  designed  to  transfer  mechanical  energy  within  the  gait  cycle  to  assist  the  person  in 
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walking,  and  c)  it  can  be  used  to  create  a  better  dynamic  biomechanical  model  (DBM)  as 
it  can  help  model  the  actions  of  the  pack-person  interface  force  requirements.  The 
objective  was  :  to  further  develop  the  Fastrak®  software  for  the  Load  Carriage  Simulator 
to  provide  phase  shift  and  amplitude  infonnation  about  the  motion  of  the  pack  relative  to 
the  mannequin. 


2.4  Software  Development 

One  objective  was  to  develop  a  protocol  that  did  not  affect  the  standardized 
approach  currently  used  in  LC  Simulator  testing.  Hence,  we  did  not  wish  to  remove  the 
Fastrak®  source  from  its  location  on  the  mannequin.  The  solution  was  to  use  a  second 
sensor  as  the  stationary  laboratory  reference  point  (0,0,0).  In  this  way,  the  reference 
(sensor  #2)  served  to  reflect  the  laboratory  coordinates  while  the  Fastrak®  source  on  the 
mannequin  measured  as  the  difference  between  laboratory  coordinates  and  the  mannequin 
movement  (source  -  sensor#2)  and  Sensor  #1  in  the  pack  measures  the  difference  between 
the  mannequin  and  the  pack’s  centre  of  gravity  (sensor#l).  A  schematic  of  this  system  is 
attached  as  Figure  2-5. 


The  relationship  of  the  source  (mannequin  motion)  to  sensor  #1  (pack  motion)  and  to 
sensor  #2  (global  reference).  The  relative  displacement  volumes  are  small  and  there  is 
no  metal  interference  in  the  acquisition  area. 

The  global  reference  sensor  is  mounted  approximately  30  cm  below  ceiling  height 
on  a  Plexiglas  rod  with  plastic  screws.  The  source  remains  on  the  mannequin’s  arm 
and  sensor  #1  remains  in  the  pack  with  the  coordinates  projected  mathematically  to 
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the  pack  centre  of  gravity.  The  protocol  for  data  collection  remains  unchanged  in 
regard  to  data  acquisition.  However,  there  is  now  additional  output  from  the  two 
sensors  into  a  data  acquisition  program.  The  main  changes  are  reflected  in  the  output 
and  analysis  programs. 

The  input  program  has  been  changed  in  the  following  ways  :  1)  it  now  allows 
multiple  trial  sets  to  be  gathered  without  reentry,  2)  it  retains  input  infonnation  from 
set  to  set  to  allow  easier  editing,  3)  it  can  collect  for  a  variable  number  of  sensors  (1  to 
4),  and  4)  the  output  file  header  has  been  revised  to  provide  additional  information  - 
number  of  sensors  and  collection  frequency  in  Hz. 

The  output  program  consists  of  two  parts:  a  graphic  representation  of  the  data  and 
an  ASCII  printout  (readable  by  Excel)  of  frequency-based  data.  The  graphs  have  been 
parceled  into  x,  y,  and  z  sheets  with  (a)  as  displacement  data  (in  mm),  (b)  as  frequency 
(in  Hz),  and  (c)  as  a  correlogram  of  phase  shift  between  the  global  and  pack 
infonnation.  Important  information  such  as  ranges,  and  principal  frequencies  are  also 
recorded.  The  peak  values  for  displacement  (a)  are  calculated  by  identifying  the  peak 
amplitudes  for  each  ten  second  interval  and  averaging  them.  The  frequency  spectrum 
(b)  for  each  coordinate  and  trial  is  calculated  from  a  fast  fourier  transform  (FFT) 
algorithm  using  5 12  samples  (approx  9.3  sec)  and  providing  amplitude  and  phase 
values  for  128  frequency  values  from  0.1  to  about  13  Hz.  The  printout  is  scaled  to 
the  highest  values  of  amplitude  within  each  trial  with  this  dominant  frequency 
reported  on  the  left  side  of  the  graph  as  the  maximal  frequency  of  the  signals  for 
global  and  pack  sensors.  The  EXCEF  printout  file  captures  both  the  FFT  and  phase 
angle  as  sine/cosine  angle  information  that  is  needed  to  return  the  signal  to  the  time 
domain.  These  data  will  be  needed  for  biomechanical  modeling  once  the  response 
characteristics  of  the  signal  are  captured  in  the  biomechanical  model  connection 
points.  The  final  data  display  (c)  is  a  correlogram  representing  a  cross-correlation 
analysis  between  the  global  signal  and  the  pack  signal.  To  create  this  graph,  the 
global  signal  is  held  constant  in  time  and  the  pack  displacement  signal  is  shifted  to  the 
right  by  one  time  step.  Then,  the  root  mean  square  values  are  calculated  and  recorded. 
This  process  is  repeated  for  a  total  of  1 .2  seconds  in  order  to  see  where  the  amplitudes 
are  largest,  a  point  that  reflects  the  phase  shift  in  seconds  between  the  two  signals.  To 
avoid  jitter  at  the  beginning  of  the  trial,  this  data  collection  routine  is  started  after  a  30 
second  delay  in  sampling.  Maximal  data  are  reported  in  terms  of  frequency  of  delay 
in  radians  and  time  of  delay  in  seconds. 

Figures  2-6,  2-7,  and  2-8  represent  the  x  (anterior/posterior),  y  (side/side)  and  z 
(up/down)  data  from  this  program.  In  the  Figures  2.6a,  2.7a  and  2.8a,  it  is  possible  to 
see  the  global  displacement  data  with  the  amplitudes  varying  from  24.06  mm,  22.48 
mm  and  47.3 1  mm  for  x,  y  and  z  displacement  respectively.  This  makes  sense  given 
the  FC  Simulator  input  function  for  gait.  The  smaller  signal  on  each  of  the  designated 
graphs  represents  the  motion  of  the  pack,  an  amount  that  is  approximately  10%  of  the 
FC  Simulator  forcing  function.  The  CTS  backpack  was  superior  to  most  packs  in  the 
database  at  1.32mm,  0.84  mm  and  0.52  mm  for  pack-person  relative  motions  in  x,  y, 
and  z  planes  (Reid  et  ah,  1999).  This  results  mean  that  the  pack  suspension  system  is 
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acting  as  a  large  damper  of  gait  motion  thus  absorbing  much  of  the  original  motion  so 
that  the  pack  sensor  is  only  receiving  a  small  percentage  of  the  original  displacement 
amplitudes.  It  can  also  be  noted  in  the  x  and  z  displacement  graphs  that  the  timing  of 
the  pack’s  peak  displacements  are  not  in-phase  with  the  mannequin’s  forcing  function. 

In  Figures  2-6b,  2-7b,  and  2-8b,  the  frequency  content  of  the  signal  is  shown.  In 
all  data  sets,  the  maximal  frequency  is  1.82  Hz.  This  is  the  frequency  of  the  forcing 
function  of  the  LC  Simulator  for  this  trial  thus  providing  a  validation  of  the  data 
calculations.  This  frequency  would  be  expected  to  dominate.  There  are  small 
harmonic  frequencies  shown  at  3.6  Hz  and  5.5  Hz,  a  factor  that  demonstrated  a 
rhythmical  nature  of  the  forcing  function.  Noise  in  the  signals  range  from  2%  to  20% 
of  signal  strength.  The  raw  displacement  data  (a)  shows  this  noise  between  data 
points  as  well.  Some  of  this  noise  can  be  removed  with  better  tuning  of  the  forcing 
function.  Another  reason  for  this  noise  was  minor  vibration  of  the  global  reference 
sensor  that  was  not  isolated  from  the  floor  during  this  sample  trial.  Lastly,  the 
Fastrak™  raw  data  signals  were  not  filtered  for  high  frequency  noise.  With 
improvements  in  tuning  and  isolation,  this  noise  will  be  reduced  to  less  than  5%.  The 
noise  does  not  interfere  with  a  clear  understanding  of  the  signal. 

In  Figures  2-6c,  2-7c,  and  2-8c,  the  correlogram  is  being  used  to  show  the  phase 
shift  or  lag  between  the  mannequin’s  motion  to  the  pack’s  motion.  In  the  x 
(forward/backward)  and  z  (up/down)  directions,  the  pack  is  out  of  phase  with  the  LC 
Simulator  by  3.02  and  3.44  radians  (173°  and  197°)  respectively.  This  is  also 
represented  by  a  delay  of  .29  seconds  and  .35  seconds  respectively.  As  the  conversion 
of  the  LC  Simulator  motion  from  1.82  Hz  to  one  cycle  period  is  0.55  seconds,  this 
represents  a  complete  out-of-phase  motion  of  the  pack.  The  y  displacement  (side  to 
side)  motion  is  almost  in  phase  with  the  LC  Simulator  function  as  its  motion  is 
maximized  at  0.66  radians  (or  37.8°)  of  phase  shift  that  can  also  be  described  as  a  time 
lag  of  .55  seconds.  Without  further  data,  it  is  not  possible  to  interpret  the  importance 
of  these  findings.  However,  it  is  speculated  that  a  particular  out-of-phase  pattern  will 
be  best  for  transfer  of  mechanical  energy  from  the  pack  to  the  person  as  speculated  by 
Kram  (2000).  We  plan  to  study  this  relationship  in  the  near  future. 
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Figure  2-6.  Relative  X  motion  between  pack  and  person:  a)  x  displacement;  b)  x  FFT 
spectrum  and  c)  x  phase  shift. 
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Relative  Amplitude  Relative  IS/btion,  (rrm) 


Figure  2-7.  Relative  Y  motion  between  pack  and  person:  a)  y  displacement;  b)  y  FFT 
spectrum  and  c)  y  phase  shift. 
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Figure  2-8.  Relative  Z  motion  between  pack  and  person:  a)  z  displacement;  b)  z  FFT 
spectrum  and  c)  z  phase  shift. 
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2.5 


Uses  of  Software 


One  immediate  use  of  the  software  is  with  the  development  of  the  new  DBM 
model.  The  resultant  motion  of  the  pack  is  required  in  order  to  model  the  constraints 
between  the  two  moving  systems.  The  software  modification  will  also  allow  us  to 
characterize  the  interface  more  accurately  so  that  more  accurate  modelling  routines  can  be 
developed. 

A  second  use  is  to  determine  if  there  are  differences  between  packs  in  terms  of  the 
mechanical  transfer  of  energy  within  the  suspension  system.  In  Kram’s  (2000) 
hypothesis,  he  suggested  that  a  natural  frequency  of  vibration  of  a  pack  is  equal  to  the 
square  root  of  k/m  where  k  is  the  spring  stiffness  and  m  is  the  mass  supported  by  the 
spring.  He  felt  that  a  pack  suspension  system  could  be  tuned  to  apply  a  steady  force  with 
moderate  force  fluctuations  and  challenged  designers  to  incorporate  compliant  suspension 
systems  into  load  carrying  systems  and  measure  the  impact  of  this  change. 

A  third  use  is  to  design  a  dynamic  suspension  system  where  the  suspension  system 
of  the  pack  can  be  tuned  to  create  a  natural  frequency  for  marching  orders.  It  would  be 
important  to  also  have  soldier  trials  to  determine  whether  there  is  a  net  savings  of 
metabolic  costs  and  whether  the  soldier  can  “feel”  that  the  pack  is  now  costing  less 
energy.  This  analysis  would  be  useful  to  detennine  the  cost-effectiveness  of  active 
suspension  systems. 
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3.0 


Construction  of  Small  Strap  Force  Gauges 


3.1  Introduction 

The  current  design  for  strap  force  sensors  is  a  dog-bone  made  out  of  a  T6061  aluminium 
base  plate  (t=1.25mm,  L=39  mm)  with  two  ninety  degree  “T”  rosette  strain  gauges 
(Measurements  Group  Inc.  M-M:  EA-13-125TG-350)  mounted  in  a  full  Wheatstone 
bridge  configuration.(Figure  3-1).  A  stress  relief  was  used  to  prevent  the  strain  gauges 
from  being  torn  from  their  cables.  Straight  pins  are  used  to  suspend  the  gauge  and  unload 
the  underlying  strap  causing  all  load  to  be  directed  through  the  gauge.  These  strap  force 
sensors  were  linear  and  robust  but  the  hook  and  pin  attachments  required  a  minimum  of  6 
cm  of  strap  length.  This  design  worked  well  for  the  shoulder  straps  and  waist  belt. 
However,  it  was  not  suitable  for  load-lifter  straps  and  hip  stabilizer  straps  because  of 
length  of  strap  available  and  the  radius  of  curvature  at  some  locations  around  the 
shoulders  and  hips. 


Figure  3-1.  Sample  “dog-bone”  strap  sensor  used  in  previous  contracts  for  shoulder 
straps  and  waist  belts. 


3.2  Purpose 

The  purpose  of  this  work  plan  was  to  construct  and  calibrate  smaller  strap  tension 
sensors  to  access  the  difficult-to-measure  straps  (i.e.,  load  lifters,  hip  stabilizers).  The 
goal  was  to  develop  eight  strap  gauges  that  could  withstand  the  same  tension  ranges  (~ 

5 ON  to  100N)  but  be  used  in  tight  locations. 

3.3  Methods 

A  series  of  aluminium  blanks  were  slit  (not  sheared)  from  1.568  mm  thick  T6061 
Aluminium  sheet.  Each  of  these  was  18  by  6.5  mm  in  size.  Mechanical  designer,  Mr. 
Gerry  Saunders  from  the  Human  Mobility  Research  Centre  (HMRC)  machined  the 
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backing  plate  shapes,  drilled  attachment  holes  at  each  end  of  the  blanks  and  cleaned  the 
metal  in  preparation  for  gauge  mounting.  Two  ninety  degree  “T”  rosette  strain  gauges 
(Measurements  Group  Inc.  M-M:  EA-13-062TT-350)  were  mounted  in  a  full  Wheatstone 
bridge  configuration  on  each  blank.  A  transparent  protective  coating  (M-Coat  A, 
Measurements  Group  Inc.)  was  used  over  the  entire  surface  to  water  proof  the  transducer. 
Physical  length  of  the  new  gauges  is  18  mm  (Figure  3-2). 


Figure  3-2.  Sample  of  the  new  miniature  strain  gauge  sewn  into  a  strap. 


Nine  miniature  strap  force  sensors  were  made.  Each  gauge  was  calibrated  with  an 
Instron™  force  measurement  device  under  static  loads  equivalent  to  the  ranges  of  force 
used  in  load  carriage  conditions.  Figure  3-3  depicts  the  testing  set-up. 


Figure  3-3.  Miniature  strap  force  sensor  calibration  in  the  Instron.  A  portable  Strain 
Indicator  was  used  to  acquire  strain  gauge  data. 
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3.4  Results  of  Calibration 


Figure  3-4  depicts  a  summary  of  the  nine  sensors  calibrated.  The  regression  R2  values 
ranged  from  0.9934  to  0.9999  with  the  correction  equations  provided  in  Table  3-1.  Only 
Gauge  #7  was  less  accurate  and  also  had  a  noisy  baseline.  The  other  gauges  are  highly 
linear  and  acceptable  for  use. 


SUMMARY  OF  SMALL  SENSORS 
CALIBRATION  CURVES  (1-9) 
(Sensor  #  M7  different) 


Microstrains 


Figure  3-4.  Summary  of  all  nine  miniature  strap  force  gauges. 


Table  3-1.  Summary  of  regression  equations  for  miniature  strain  gauges. 


Gage 

Regression  Equation 

R2 

#1 

y=  6.1355x  +16.761 

0.999 

#2 

y=  6. 1895x+  1.1204 

0.9999 

#3 

y=5.9705x- 6.2452 

0.9997 

#4 

y=  6.0831x  +  3.1883 

0.9999 

#5 

y=6.1236x- 0.8172 

0.9995 

#6 

y=  6.015  x-  4.8406 

0.9998 

#7 

y=6.7419x  + 37.918 

0.9934 

#8 

y=6.2446x-  1.0136 

0.9999 

#9 

y=5.9544x- 4.3228 

0.9999 
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3.5 


Discussion  and  Conclusions 


The  miniature  strap  gauges  will  be  suitable  for  use  in  most  locations.  It  must  be  noted 
however,  that  these  gauges  have  not  satisfied  three  requirements:  they  cannot  be  mounted 
quickly  for  field  trials,  they  require  the  strap  to  be  floating  free  of  underlying  material  (a 
difficult  task  given  the  nature  of  the  pack  design),  and  they  may  be  more  susceptible  to 
damage  from  bending  of  the  sensor.  These  design  problems  will  be  corrected  with  the 
new  strap  sensor  design  that  is  part  of  the  Phase  II  contract.  In  addition,  in  the  new  sensor 
design,  there  would  be  no  need  to  modify  the  strap  (cut  into  or  sew  onto)  and  it  is  likely 
that  a  sensor  could  be  captured  directly  by  a  portable  data  acquisition  system  with  no 
requirement  for  an  amplifier  (saving  power/weight  and  reducing  cost). 

In  conclusion,  and  despite  these  possible  limitations  mentioned  above,  these  nine  strap 
sensors  are  only  1 8mm  in  length,  highly  linear,  robust  and  a  substantive  improvement  in 
size  from  the  previous  “dog-bone”  sensor  design. 

4.0.  Trifilar  Pendulum:  Experimental  Determination  of  the  Moment  of 
Inertia 

4.1  Introduction 

One  requirement  of  dynamic  biomechanical  modeling  (DBM)  is  to  know  the  location  of 
the  centre  of  gravity  and  the  moment  of  inertia  of  a  pack.  It  is  possible  that  two  packs  with 
the  same  load  and  a  completely  different  geometry  can  have  the  same  centre  of  gravity 
and  moment  of  inertia.  Due  to  the  non-uniform  nature  of  a  backpack,  pinpointing  the 
centre  of  gravity  can  be  difficult.  While  the  packs  ideally  would  be  loaded  uniformly,  it  is 
generally  not  possible.  Attempts  can  be  made  to  centre  the  load  or  have  it  in  the  same 
position  relative  to  the  user  but  the  non-rigid,  non-uniform  shape  of  the  compartments 
make  it  impossible. 

Mathematically  calculating  the  moment  of  inertia  of  a  non-unifonnly  shaped  and 
weighted  object  requires  complex  equations.  With  the  proper  tools,  this  process  is 
simplified  by  physically  calculating  the  centre  of  gravity  and  moment  of  inertia.  The  most 
logical  method  is  to  construct  a  trifilar  pendulum  to  make  these  calculations. 

The  trifilar  pendulum  provides  a  relatively  simple  and  quick  method  of  locating  the  centre 
of  gravity  and  the  moment  of  inertia  of  non-homogenous  objects.  Its  applications  go  far 
beyond  those  of  the  backpack  project.  It  can  be  used  to  measure  properties  of  objects 
other  than  backpacks  that  are  worn  or  carried  by  soldiers.  In  addition,  the  pendulum  is 
large  and  stable  enough  to  measure  centre  of  gravity  and  moment  of  inertia  of  the  human 
body  in  any  position. 

4.2  Theoretical  Basis 
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According  to  Webster’s  Dictionary,  moment  of  inertia  is  ‘the  ratio  of  the  torque  applied  to 
a  rigid  body  free  to  rotate  about  a  given  axis  to  the  angular  acceleration  thus  produced, 
about  that  axis  and  equal  to  the  sum  of  the  products  of  each  element  of  each  mass  by  the 
square  of  its  distance  from  the  given  axis.  Moment  of  inertia  is  to  rotation  as  mass  is  to 
linear  movement.  The  equation  for  this  is: 


T=Ia 


The  word  trifilar  consists  of  ‘tri’,  or  three  and  ‘filar’,  or  filaments;  a  pendulum  that  is 
hung  with  three  filaments.  It  resembles  a  tire  swing  with  the  exception  of  the  insertion 
points  of  the  top  of  the  filaments.  Instead  of  converging  at  one  point,  each  filament 
remains  completely  vertical  with  the  insertion  point  directly  above  the  point  of  insertion 
on  the  pendulum  plate. 


The  trifilar  pendulum  consists  of  a  round 
plate  of  a  homogenous  material,  suspended 
in  the  horizontal  plane  by  three  cables  of 
equal  length  equidistantly  placed  from  the 
centre  of  the  plate.  When  the  pendulum  is 
twisted  around  a  vertical  axis  running 
through  the  centre  of  the  plate,  it  oscillates 
with  a  period  dependent  on  the  following 
factors:  the  length  of  the  cables,  the  radius  of 
the  plate,  the  mass  of  the  plate,  the  local 
force  of  gravity,  and  the  moment  of  inertia  of 
the  plate  [See  equation  1,  Appendix  A]. 


i  .  ■  ,  ,  ,,  ,  ,  ■  ,  Figure  4-1.  Trifilar  Pendulum  with  Pack 

When  an  object  is  placed  on  the  plate,  with  ° 

its  centre  of  gravity  (c  of  g)  above  the  c  of  g 

of  the  plate,  the  resultant  object/plate  combination  can  be  treated  as  a  new  pendulum. 
Since  all  factors  other  than  mass  remain  constant,  the  change  in  the  period  of  the  new 
pendulum  is  related  to  differences  in  the  mass  of  moments  of  inertia  of  the  old  and  new 
pendulums.  The  mass  can  be  easily  measured  and  its  effect  calculated,  and  the  moment  of 
inertia  of  the  object  can  be  detennined  using  a  mathematical  formula  [See  equations  2  and 
3,  Appendix  A]. 


The  design  is  based  on  two  smaller  trifilar  pendulums  that  are  presently  used  for  2nd  and 
3ld  year  engineering  labs  at  Queen’s  University.  Although  these  pendulums  are  too  small 
for  measuring  backpacks,  they  represent  a  proven  design  that  can  be  scaled  up  for  larger 
objects.  Several  modifications  were  made  to  increase  the  stability  and  safety  for 
measuring  properties  of  humans. 
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4.3  Strengths  and  Limitations  of  Design 

There  are  several  strengths  of  the  use  of  the  tri  filar  pendulum.  It  requires  very  simple 
calculations  that  can  be  done  on  a  basic  calculator.  Once  the  set-up  is  installed,  it  provides 
quick  results  and  is  easy  to  use.  It  is  relatively  inexpensive  and  upkeep  and  maintenance  is 
minimal.  The  equipment  is  durable  and  should  stand  up  to  years  of  use.  It  makes  an 
excellent  teaching  tool  as  it  has  a  wide  variety  of  applications  and  is  quite  versatile.  It  is 
more  accurate  than  many  of  the  other  commonly  used  methods. 

There  are  several  limitations  to  the  design  as  well.  The  trifilar  pendulum  has  been 
optimally  designed  for  objects  no  larger  than  0.89m  wide,  although  as  long  as  the  entire 
mass  of  the  object  is  supported  by  the  plate  any  length  is  acceptable.  Due  to  the  geometry 
of  the  pendulum,  the  accuracy  increases  with  weight.  Therefore  objects  lighter  than  13.61b 
have  +/- 1 0%  error.  Human  error  of  the  scale  reading  is  variable.  Of  the  mechanical  errors 
involved,  one  must  consider  the  inherent  errors  in  the  scale  readout  (+/-0.45kg), 
dimensional  tolerances  of  the  plate  and  irregularities  of  the  plate  lamination,  accuracy  of 
the  filament  connections,  accuracy  of  the  placement  of  the  scales,  horizontal  movement  of 
the  plate  during  rotation,  and  shifting  of  the  load  during  testing.  Non-mechanical  errors 
attributed  to  the  operators  include  inaccurate  scale  reading,  errors  in  timing  and  complete 
rotation  detennination,  inadvertent  horizontal  movement  and  the  positioning  for  the  angle 
of  rotation  of  the  plate. 

4.4  Design  Requirements 

Several  design  requirements  had  to  be  considered,  they  are  as  follows: 

1 .  Accurate  and  precise  results 

The  current  methods  used  for  finding  the  moments  of  inertia  of  packs  require  a  great  deal 
of  guesswork  and  complex  calculations.  The  trifilar  pendulum  design  provides  more 
mathematically  accurate  results  and  has  a  high  level  of  repeatability. 

2.  Location  (requires  a  high  ceiling  for  mounting) 

Due  to  the  long  length  of  the  filaments,  installation  location  was  important.  Shortening  of 
the  filaments  due  to  height  restraints  would  change  the  properties  of  the  pendulum.  The 
importance  of  a  non-permanent  assembly  was  also  important,  as  the  apparatus  is  located 
in  a  classroom.  The  design  allowed  for  the  top  cable  attachments  to  be  pennanently  fixed 
to  the  ceiling  while  the  pendulum  plate  and  lower  plate  could  be  easily  moved.  The  cables 
can  then  be  hooked  up  onto  the  ceiling  out  of  the  way  as  the  rest  of  the  assembly  is  stored 
in  a  closet  until  it  is  needed. 

3.  Wide  range  of  applications/different  objects 
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The  pendulum  was  to  have  more  applications  than  just  use  with  backpacks.  The  size 
should  be  optimal  for  use  with  a  variety  of  objects.  The  size  required  for  packs  was  used 
as  the  main  factor,  but  consideration  was  given  for  larger  and  slightly  smaller  objects  as 
well.  The  pendulum  can  be  scaled  down  if  needed  for  smaller  objects. 

The  maximum  weight  allowable  for  objects  is  163.6kg. 

4.  Ease  of  use 

Ease  of  use  was  a  key  issue  as  the  Queen’s  pendulum  is  mounted  in  a  classroom.  Ease  of 
use  also  encompasses  ease  of  maintenance  and  durability.  Modifications  to  the  initial 
design  were  made  to  make  the  attachment  of  the  plate  to  the  filaments  easier.  With 
heavier  objects  resting  on  the  plate,  it  was  difficult  to  lift  the  plate  to  attach  the  filaments 
after  finding  the  centre  of  gravity  of  the  object.  A  method  to  raise  and  lower  the  plate 
without  unhooking  the  filaments  was  designed.  The  pendulum  can  be  easily  used  by  one 
or  two  people. 


5.  Relative  low  cost 

Relative  low  cost  was  important,  though  not  a  limiting  factor.  The  simplicity  of  the 
pendulum  theory  allowed  for  this  consideration  to  be  met  easily. 

6.  Safety 

The  safety  of  use  is  an  important  factor  as  the  pendulum  was  to  be  used  for  finding 
inertial  properties  of  humans  as  well  as  inanimate  objects.  Several  modifications  were 
made  from  the  original  design  to  increase  the  degree  of  safety.  The  addition  of  a  spindle 
bearing  assembly  in  the  cen  tre  of  the  pendulum  plate  limited  the  plate ’s  horizontal 
movement  and  assisted  with  smooth  rotation.  Error  effects  of  this  spindle  assembly  have 
not  yet  been  determined. 


4.5  Final  Design  Parts  List 


Details  of  the  final  design  are  provided  in  the  Parts  List  (Table  4-1).  The 
explanation  that  follows  will  comprise  the  final  design. 
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Table  4-1  -  Parts  List 


Itfffrnirialif 

Name 

Description 

Material 

1 

Ceiling  Template 

Square  or  circle  approx,  dia/width  40" 

Plywood  or  similar 
3/4" 

4 

Ceiling  Template  Fasteners 

According  to  ceiling  material 

Steel 

1 

Rotating  Plate 

Circle  dia.  40",  6  holes  drilled  equidistant  -  dia.  35",  mass=261b 

Fir  GIS  Plywood  3/4" 

1 

Lower  Plate 

Circle  dia.  45" 

Fir  GIS  Plywood  3/4" 

3 

3/16"  Steel  Cable 

3/32"  dia.,  135"  effective  length 

Stainless  Steel  Cable 

6 

Wire  Crimpers 

3/32"  cable  size 

Steel/Aluminum 

6 

Clevis 

Plain  Stem  Yoke  End  Clevis  dia  1/2",  length  2-1/2" 

Steel  (painted) 

6 

Clevis  Pin 

1/2"  dia,  1-1/4"  effective  length 

Stainless  Steel 

6 

Clevis  Pin  Bushing 

1/2"  I.D,  3/4"  O.D.,  9/16"  length 

Brass 

6 

Cotter  Pin 

5/32"  dia,  1-1/2"  length 

Stainless  Steel 

9 

Washer 

1/2"  Screw  Size  (9/16"  I.D.,  1-3/8"  O.D.) 

Zinc  Plated  Steel 

1 

Centre  Post 

Total  length  3”,  2”  dia  flange,  !4”  thick,  machined  to  0.6693"  dia,  2- 
3/4"  length, 

Steel 

2 

Fasteners  for  Centre  Post 

Woodscrew,  % ”  length 

Steel 

1 

Centre  Socket 

1.6535"  I.D.,  1.9035  O.D.,  2-1/4"  height 

Steel 

1 

Center  Socket  Flange 

3-1/2"  dia.,  1/4"  thickness 

Steel 

1 

Bearing 

Needle  roller  bearing  with  flanges,  without  inner  ring  SKF 

RNAO  17x25x1 3 

StSteel 

1 

Bearing 

Maintenance-Free  Spherical  Plain  Bearing  SKF  GE25C 

StSteel 

5 

Fasteners  for  Socket 

Woodscrew,  %”  length 

Steel 

2 

Pegs/Rollers 

1/2"  dia,  3"  length 

Steel 

3 

Head  Swivel  Screws 

9/32"  dia,  2"  length,  Techni-Tool  #81-10 

Steel 

3 

Threaded  Inserts 

1/2"  I.D.,  3/4"  length 

Steel 

3 

Bathroom  Scales 

Danze  63-8083-4  (+/- 1  lb) 

Laser  Centering  Device 
(optional) 
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4.6  Validation  Testing  to  Find  Centre  of  Gravity 

For  the  trifilar  pendulum  to  yield  the  correct  moment  of  inertia  for  an  object  being  tested 
on  the  pendulum,  it  is  necessary  that  the  c  of  g  of  the  object  be  directly  above  the  c  of  g  of 
the  plate,  so  the  centroid  of  the  object  and  the  plate  are  in  line.  If  the  object  is  regular,  it  is 
easy  to  calculate  the  c  of  g  mathematically.  However,  if  the  object  is  not  symmetrically 
shaped,  or  has  an  uneven  mass  distribution,  it  can  be  quite  difficult  to  locate  the  c  of  g 
mathematically. 

Therefore,  a  system  was  required  that  would: 

1.  Allow  the  centre  of  gravity  of  an  object  to  be  determined  easily 

2.  Aid  the  operator  in  centering  the  object  on  the  plate. 

4.6.1  Proposed  design  being  tested 

The  proposed  design  consisted  of  a  plywood  plate,  0.5  lm  in  radius,  with  three  drilled 
holes  at  0.44m  from  the  centre.  The  holes  formed  the  three  vertices  of  an  equilateral 
triangle.  Weigh  scales  were  used  for  measuring  the  forces  at  the  three  drilled  points. 

These  scales  were  calibrated  using  known  weights  and  found  to  have  a  precision  of  +/- 
0.45kg.  Force  was  transmitted  using  a  ball  bearing  to  ensure  that  the  force  was  being 
transmitted  at  a  single  point.  From  statics,  it  is  possible  to  derive  the  relationship  between 
the  c  of  g  of  the  object  and  the  vertical  force  on  the  three  locations  by  the  summation  of 
moments. 

4.6.2  Methodology  of  Testing 

The  system  was  calibrated  at  fifteen  distinct  locations  on  the  plate.  Following  the  theory 
underlying  the  device  (that  the  moment  of  the  c  of  g  is  balanced  by  the  moments  created 
by  the  three  forces),  only  locations  (points)  inside  an  imaginary  triangle  formed  by  the 
three  points  were  used.  If  locations  outside  the  triangle  were  used,  the  plate  would  tip,  and 
the  results  would  be  useless. 

The  weigh  scales  were  zeroed  to  remove  the  effect  of  the  weight  of  the  plate.  A  13.6kg 
weight  with  a  known  centre  of  gravity  was  then  placed  such  that  its  centre  of  gravity 
exactly  coincided  with  one  of  the  fifteen  calibration  locations.  The  scales  were  then  read, 
and  the  results  recorded.  The  object  was  then  removed,  the  scales  re -zeroed,  and  the 
object  was  placed  on  the  second  calibration  location.  This  process  was  repeated  until  all 
fifteen  calibration  locations  had  been  tested.  The  entire  trial  was  then  repeated  two  more 
times. 

The  raw  data  (location,  measured  weights)  were  then  analyzed,  using  the  predictive 
equations  for  the  X  and  Y. 

4.6.3  Statistical  analysis  of  results 

If  the  c  of  g  finder  performed  perfectly,  there  would  be  a  1:1  correlation  between  the 
measured  X  and  Y  values,  and  the  calculated  X  and  Y  values.  There  was  a  high 
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correlation  between  the  actual  position  of  the  object  and  the  calculated  c  of  g  from  the 
testing.  On  the  X  dimension,  linear  regression  yielded  a  relationship  of  0.98: 1  with  an  r 
squared  value  of  >0.99.  On  the  Y  dimension,  linear  regression  yielded  a  relationship  of 
0.99: 1  with  a  r  squared  value  of  >0.99.  This  near  linear  relationship  confirms  that  the 
theory  relating  the  weights  at  the  three  triangle  vertices  to  the  location  of  the  object  is 
correct.  This  relationship  can  be  observed  in  Appendix  D. 

4.6.4  Error  Analysis 

The  error  ranged  from  -3.3  to  2.3  cm  on  the  X  dimension,  and  from  -4.8  to  2.8  cm  on  the 
Y  dimension  (fig  8,  9).  The  magnitude  of  the  error  does  not  appear  to  be  affected  by  the 
distance  from  centre  (fig  10,  11).  The  frequency  of  the  magnitude  of  the  errors  appears  to 
follow  a  normal  distribution  pattern  (fig  8,9). 

The  total  error  (ET)  of  the  system  is  calculated  with  the  Pythagorean  formula.  When  each 
trial  was  completed  separately,  the  average  ETfor  the  three  trials  was  1.73,  1.4,  and  1.5 
cm.  Since  the  errors  appear  to  cluster  around  zero,  (fig  10,  1 1)  it  was  hypothesized  that 
the  error  could  be  reduced  by  averaging  the  results. 

For  each  location,  the  weights  at  the  three  vertices  were  measured  three  times,  once  per 
trial.  For  each  location,  therefore,  there  were  three  calculated  X  and  three  calculated  Y 
values.  The  X  values  were  averaged  to  yield  a  new  (average)  X  value,  and  the  Y  values 
were  averaged  to  yield  a  new  (average)  Y  value. 

The  average  ETof  the  new  X,  Y  pairs  was  only  0.36  inches  (Appendix  E). 

This  result  can  also  be  seen  on  the  position  vs.  error  charts,  where  the  average  values  tend 
to  be  more  accurate  (less  error)  than  the  non-averaged  values. 

The  level  of  accuracy  of  this  device  could  be  of  concern  for  some  applications.  The 
average  absolute  total  error  is  1.55cm  (3.5%  of  the  radius),  and  the  maximum  absolute 
error  was  4.83cm  (10%  of  the  radius).  It  appears  that  the  accuracy  can  be  increased 
significantly  by  averaging  the  results  of  multiple  trials.  Indeed,  averaging  three  trials 
yields  an  averaged  absolute  error  of  0.91cm  (down  from  1.55cm)  which  is  2.5%  of  the 
radius.  With  this  method,  the  maximum  absolute  error  drops  to  1.78cm  (down  from 
4.83cm)  which  is  4%  of  the  radius  of  the  plate. 

The  accuracy  of  the  device  should  increase  as  the  weight  of  the  object  being  measured 
increases  due  to  the  precision  of  the  scales.  The  error  is  +/-  0.45kg,  therefore  for  a  13.6kg 
object  the  error  is  +/-  3.3%,  and  for  a  45.5kg  object  the  error  is  +/-!%. 


4.7  Validation  Testing  for  the  Moment  of  Inertia 

In  order  to  validate  the  accuracy  and  precision  of  the  pendulum,  a  known  mass  with 
homogenous  weight  distribution  and  an  easily  calculated  center  of  mass  and  moment  of 
inertia  was  chosen  as  the  object.  The  object  was  rectangular  solid  with  a  mass  of  17.1  kg 
The  object  was  tested  using  a  modification  of  the  “Backpack  Protocol”  in  the  following 
section,  although  only  one  axis  was  tested  due  to  the  shape  of  the  object.  Instead  of  the 
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usual  three  repetitions  per  axis,  ten  were  perfonned  on  the  single  axis  to  determine 
whether  or  not  there  were  order  effects.  After  ten  repetitions,  the  object  was  removed  and 
the  apparatus  was  reset.  The  test  was  repeated  for  another  ten  repetitions,  reset,  and 
repeated,  for  a  total  of  three  sets  of  data. 

An  ANOVA  two-factor  without  replication  analysis  was  run.  The  p-value  between  the 
three  sets  was  0.83  and  the  p-value  between  ordered  repetitions  was  0.35.  A  p-values  that 
is  greater  than  0.05  is  considered  to  show  no  effect  of  order  of  repetitions  or  sets.  No 
effects  were  shown. 

The  mean  measured  time  for  10  rotations  was  19.863  s  with  standard  deviation  of  0. 1 13. 
The  standard  error  is  0.02  (Appendix  D). 


4.7.1  Initial  Set-up 

A  suitable  room  was  chosen  to  accommodate  the  pendulum.  The  length  of  the  cables  was 
calculated  precisely  and  was  not  changed,  as  this  would  have  effected  the  mathematical 
calculations.  The  length  of  cables  was  the  optimal  value  for  the  size  and  weight  of  the 
suspended  plate. 

The  ceiling  template  was  affixed  to  the  ceiling  with  proper  fasteners  that  exceed  the  load 
that  the  pendulum  is  rated  to  hold.  Care  was  taken  to  ensure  that  the  plate  was  level. 
Three  cables  were  attached  to  the  ceiling  plate  with  the  cotter  pins. 
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4.7.2  Test  Protocol 


1 .  Place  the  bottom  plate  on  the  floor  below  the  ceiling  plate.  Arrange  the  three  scales  on 
the  plate  with  the  dials  facing  out  and  their  rubber  feet  in  the  corresponding  indents  on 
the  plate. 

2.  Clean  the  shaft  on  the  underside  of  the  pendulum  plate  to  ensure  it  is  free  of  dust  and 
debris.  Insert  it  into  the  socket  on  the  bottom  plate.  This  must  be  done  carefully  as  the 
socket  contains  a  bearing. 

3.  Tare  the  scales. 

4.  Attach  the  cables  to  the  turnbuckles  on  the  pendulum  plate.  Reposition  the  bottom 
plate  in  order  to  line  up  the  swivel  head  screws  with  their  respective  indents  on  the 
tops  of  the  scales. 

5.  Adjust  the  swivel  head  screws  to  a  unifonn  height,  ensuring  that  the  full  weight  of  the 
plate  is  on  them.  The  cables  should  be  slack. 

6.  Tare  the  scales. 

7.  Place  the  load  on  the  pendulum  plate.  Position  it  such  that  each  scale  is  weighted 
evenly. 

8.  Gently  lower  the  swivel  head  screws  until  the  cables  are  tensioned  and  the  plate  is 
suspended.  If  the  load  shifts,  lower  the  plate  and  repeat  the  centering  steps. 

9.  Insert  the  angle  measurement  pegs. 

10.  Rotate  the  pendulum  plate  slowly  until  it  is  displaced  ten  degrees. 

11.  Release  the  pendulum  plate  in  conjunction  with  the  beginning  of  timing. 

12.  Record  the  time  required  to  complete  ten  full  rotations. 

13.  Repeat  steps  10  through  12  three  times  and  average  the  final  result. 

14.  Take  a  digital  photo  from  as  close  to  directly  above  the  pendulum  plate  as  possible. 

15.  Raise  the  swivel  head  screws  until  the  plate  is  evenly  supported  on  the  scales. 

16.  Reposition  the  object  for  calculations  about  a  different  axis. 

17.  Repeat  steps  7  through  14. 

18.  Reposition  the  object  for  calculation  about  the  final  axis,  steps  7  through  14. 
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4.7.3  Calculations 


The  governing  equation  used  to  calculate  the  moment  of  inertia  about  each  axis  is: 


Where  : 

/  =  mass  moment  of  inertia  (kg  ■  m2) 
m  =  mass  of  plate/object  combination  (kg) 

T  =  period  of  plate/object  combination  (5) 
mp  -  mass  of  plate  (kg) 

Tp  =  period  of  plate  (5) 
r  =  effective  radius  of  plate  (m) 
l  =  length  of  filaments  (m) 

Details  can  be  found  in  Appendices. 

4.7.4  Location  of  the  Centre  of  Gravity  of  a  Backpack 


Figure  4-2.  To  locate  the  centre  of 
gravity,  the  photos  taken  after  each 
pack  position  are  manipulated  as 
seen.  The  point  at  which  the 
bisecting  lines  cross  in  the  centre  is 
the  position  of  the  centre  of  gravity 
of  the  pack  in  that  plane. 


4.7.5  Conclusions 


The  trililar  pendulum  design  has  proven  to  be  a  simple  and  user  friendly  apparatus  to 
determine  the  centre  of  gravity  and  moment  of  inertia  of  irregularly  shaped  objects.  Most 
of  the  design  considerations  have  been  met,  although  the  accuracy  of  the  design  will  be 
improved  with  testing  of  the  new  model,  outlined  in  the  detailed  drawings. 
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5.0  Automation  of  the  Load  Carriage  Compliance  Tester 


5.1  Introduction 

Immobilisation  of  the  torso  has  been  shown  to  increase  a  person’s  energy  expenditure  in 
walking  at  moderate  speeds  by  about  10%.’  Energy  cost  was  more  than  twice  that  of 
carrying  a  5kg  weight  about  the  waist.  When  applied  to  load  carriage,  the  upshot  is:  the 
greater  the  resistance  to  nonnal  body  movement,  the  greater  the  baseline  energy  cost  for  a 
load  carriage  system.  As  well,  good  torsional  compliance  is  necessary  to  assure  free 
relative  motion  between  the  shoulders  and  hips  during  agility  activities.  The  rigidity  of  a 
LCS  is  inversely  correlated,  (r">0.86)  to  a  number  of  performance  parameters  such  as: 
comfort  of  the  user,  a  wearer’s  ability  to  perform  whole  body  motions,  and  the  ability  to 
perfonn  tasks  requiring  upper  arm/shoulder  mobility  such  as  moving  ones’  arms  to  the 
front  or  overhead. 

The  compliance  tester  was  originally  developed  to  assess  the  stiffness  of  equipment  worn 
on  the  torso  using  quasi-static  motion.  Movement  of  the  test  torso  was  achieved  using  a 
series  of  hand  operated  pulleys  and  cables.  This  required  an  operator  to  attempt  to 
produce  a  constant  rate  of  displacement  which  introduces  some  variability  in  the 
methodology  and  confined  testing  to  quasi  static  conditions. 

Conversion  to  automated  programmable  motion  control  creates  a  number  of  functions  not 
possible  with  the  previous  system.  These  include: 

•  Highly  repeatable  motion  profiles  independent  of  the  operator 

•  Determination  of  the  system  stiffness  under  dynamic  conditions 

•  Provide  quantitative  validation  of  models  of  load  carriage  devices 

•  Determine  the  frequency  response  of  LC  suspension  systems 

•  Create  a  automated  test  cell  requiring  minimal  operator  expertise  and  low  cost 
for  possible  sale  to  support  other  countries  modelling  efforts 


5.2  Purpose 

The  purpose  for  this  design  modification  is  to  support  the  development  of  a 
biomechanical  model  of  human  load  carriage.  The  device  will  initially  furnish  input 
values  for  the  dynamic  stiffness  variables  during  initial  modelling  stages.  At  later  stages 
of  the  model  development,  the  apparatus  can  be  used  to  validate  the  model  by  comparing 
the  model  predicted  stiffness  responses  to  a  different  range  of  motion  profiles. 
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5.3  Methodology 


5.3.1  Hardware  Modifications 

The  test  apparatus  consists  of  a  human  torso  shape  with  independent  thorax  and  pelvic 
sections.  These  sections  rotate  and  pivot  relative  to  each  other  through  the  required 
ranges  of  motion.  Previously,  this  motion  was  achieved  by  using  a  combination  of  a  hinge 
with  an  oil  lite  bearing  to  allow  flexion,  and  a  thrust  bearing  to  allow  relative  twist.  The 
hinge  joint  required  the  operator  to  loosen  a  number  of  locking  nuts  on  two  collars  and 
rotate  the  hinge  within  the  torso  and  then  retighten  the  locking  collars.  This  hinge  had  to 
be  aligned  parallel  to  the  path  of  motion  of  the  overhead  slider  that  applies  the  force  to  the 
top  of  the  torso.  Even  slight  misalignments  cause  increased  friction  in  the  joint.  This 
effect  had  to  be  compensated  for  by  recording  the  unloaded  force/displacement  behaviour 
of  the  test  apparatus  and  subtracting  this  effect  from  the  results  of  the  loaded  test.  These 
bearings  have  been  replaced  with  a  custom  spherical  bearing  that  allows  both  twist  and 
flexion  and  will  not  require  multiple  realignments  by  the  operator. 

Prior  to  this  upgrade,  motion  was  generated  manually  using  cables  and  pulleys. 
Subsequent  to  this  upgrade,  the  motions  of  forward  and  lateral  flexion  are  now  fully 
automated.  Linear  cyclic  motion  is  generated  using  a  programmable  motion  control 
system.  This  consists  of  a  linear  actuator  with  a  1.25m  travel:  ERB50-BLTRA20- 
BSR1250-A,  a  step  and  direction  servo  drive:  TQ10SD,  a  servo  motor:  SE233BE- 
KFLG10  and  a  PC  controller  card:  OEM-AT6400.  All  components  are  products  of 
Parker  Hannifin  Corporation,  Automation  Division,  Cleveland,  Ohio.  The  resistance  of  a 
LC  system  to  torsion,  forward  flexion,  and  medial/lateral  flexion  is  measured  in  three 
separate  tests.  Originally,  moments  applied  during  forward  and  lateral  flexion  testing 
were  measured  by  calculating  the  difference  between  the  tension  in  two  cables.  This 
required  two  single  axis  load  cells.  This  has  been  replaced  with  a  single  bending  beam 
load  cell  which  measures  the  applied  load  directly. 

Figures  5-1  illustrates  the  redesigned  compliance  tester  configured  for  lateral  flexion 
testing.  The  range  of  motion  required  for  the  lateral  flexion  test  is  shown  in  Figure  5-2. 

The  long  travel  of  the  linear  actuator  can  be  used  to  create  a  forward  flexion  of  up  to  45°. 
This  range  is  consistent  with  the  forward  flexion  angle  used  in  previous  quasi-static 
testing.  Figure  5-3  shows  the  configuration  and  range  of  displacement  of  the  test 
apparatus  for  forward  flexion  testing.  In  addition  to  full  range  forward  flexion  testing,  it  is 
possible  to  reproduce  the  small  range  forward  flexion  pattern  typical  of  walking  and  cycle 
the  test  for  a  fixed  duration  of  time.  This  gait  pattern  of  dynamic  forward  flexion  could 
be  run  multiple  cycles  and  the  forward  flexion  work  done  by  a  user  could  be  calculated. 
Repeating  this  type  of  cyclic  test  for  lateral  flexion  would  permit  calculation  of  the  lateral 
flexion  work  done.  If  a  cyclic  torsional  test  was  also  perfonned,  the  total  energy  cost  of 
the  load  carriage  system  to  the  user  could  be  estimated.  Although  originally  intended  to 
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test  rucksacks,  webbing  and  load  carriage  vest  combinations,  the  design  of  the 
Compliance  tester  allows  many  different  types  of  equipment  to  be  tested  in  this  way, 


Figure  5-1.  Compliance  Tester  Configured  for  Dynamic  Lateral  Flexion  Testing: 
a)  Top  View  and  b)  Front  View. 
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Figure  5-2.  Range  of  Motion  during  Lateral  Flexion  Testing 
(Default  frequency,  1.8  Hz). 


including  different  combinations  of  clothing  layers  and  or  body  carried  equipment  can  be 
assessed.  Physiological  measures  can  capture  the  resultant  whole  body  costs  of  a  “soldier 
system”  while  dynamic  compliance  testing  based  on  human  motions  can  help  identify 
where  that  work  is  being  done. 
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Actuator  range  of  Motion:  1220mm 
Forward  Flexion  0  -  950mm  (45°) 
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Swivel  Joint  -  on 
bending  load  cell  holds 
linear  bearing. 


Spherical  Joint  - 
allows  bend  and  twist 


Figure  5-3.  Compliance  Tester  Configuration  and  Available  Range  of  Motion  for  forward 
Flexion  Testing 
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5.2.2  Software  Modifications 

Data  acquisition,  motion  control  and  the  user  interface  were  realized  with  LAB  VIEW 
V6.1  (National  Instruments®)  allowing  integration  of  these  three  functions.  This 
integration  greatly  reduces  the  number  of  operator  interventions  required  during  testing. 
Motions  for  forward  and  medial/lateral  flexion  have  been  fully  automated  at  this  time.  An 
appropriate  linkage  to  automate  trunk  twist  has  not  been  designed  as  yet,  although  the 
data  acquisition  and  user  interface  for  torsional  testing  have  been  upgraded.  A  typical 
user  screen  for  data  acquisition  and  display  during  testing  appears  as  Figure  5-4. 
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Figure  5-4  Operator  Screen  (Graphical  User  Interface) 
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5.3  Conclusions  and  Recommendations 


The  compliance  tester  has  been  successfully  upgraded  to  perform  programmable  dynamic 
compliance  testing  in  support  of  ongoing  biomechanical  modeling  of  human  load 
carriage.  Both  the  lateral  and  the  forward  flexion  motions  have  been  successfully 
automated  with  a  corresponding  upgrade  to  the  user  interface  that  permits: 

•  Highly  repeatable  motion  profiles  independent  of  the  operator 

•  Determination  of  lateral  and  forward  stiffness  under  dynamic  conditions 

•  Quantitative  validation  of  analytical  models  of  load  carriage  devices 

•  Allow  detennination  of  the  frequency  response  of  LC  suspension  systems 

•  Maintained  backwards  compatibility,  allowing  direct  comparison  to  measures 
taken  with  the  quasi-static  design.  This  maintains  the  integrity  of  the  existing 
database. 

Additionally,  the  replacement  of  the  thrust  and  hinge  bearings  with  a  single  spherical  joint 
that  pennits  both  twist  and  flexion,  reduces  the  possibility  of  poor  alignment.  This 
change  also  greatly  reduces  the  level  of  required  operator  intervention  and  expertise. 
Much  of  the  objective  of  creating  an  automated  test  cell  that  requires  minimal  operator 
expertise  has  been  met. 

The  final  modification  required  is  the  design  and  addition  of  components  that  will  permit 
the  use  of  the  existing  linear  actuator  to  create  the  rotational  motion  required  for  torsional 
stiffness  testing.  This  will  create  a  complete  low  cost  test  apparatus  for  possible  sale  to 
support  other  countries  modelling  efforts. 

Finally,  physiological  measures  can  capture  the  whole  body  energy  costs  of  a  “soldier 
system”  while  the  dynamic  compliance  testing  based  on  human  motions  can  now  begin  to 
identify  where  that  work  is  being  done. 


6.0  Creation  of  a  Full  Body  Mapping  of  the  Mannequin 

6.1  Introduction 

The  Phase  I  static  biomechanical  model  used  only  pack  geometry  (the  shoulder  shape, 
strap  angles,  tension  and  friction  ratio  of  T1/T2  straps)  to  calculate  the  shoulder  reaction 
forces  (Stevenson  et  al.,  1995).  The  lumbar  reaction  force  was  based  on  the  lean  angle 
and  the  resultant  force  needed  to  maintain  static  equilibrium.  A  static  model  of  a  waist 
belt  was  added  in  Phase  II.  This  belt  was  modeled  as  a  partial  half  cone  where  the  slope  of 
the  cone  represented  the  anatomical  slope  of  the  hips.  The  hoop  stress  equation  was  used 
to  detennine  the  net  compressive  force  generated  by  the  tension  in  the  waist  belt.  Net 
compressive  force  could  then  be  divided  into  the  lifting  and  the  frictional  components. 
Alan  Rigby,  in  his  M.Sc.  thesis,  attempted  to  include  additional  pack  elements  (i.e.  load 
lifter  straps)  to  improve  the  model  and  validate  measured  waist  reaction  forces  with  a  load 
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cell.  He  found  widely  ranging  frictional  forces  and  variable  compliance  of  the  suspension 
system  based  on  strap  tensions.  His  conclusion  was  that  properly  modeling  the  pack- 
person  interface  forces  was  essential.  This  thesis  served  to  point  the  way  to  a  new 
dynamic  modeling  approach  that  deals  with  the  pack-person  interface  of  dynamic 
suspension  systems.  It  also  pointed  to  the  need  to  find  better  ways  to  detennine  the 
nature  of  the  pack-person  interface. 

6.2  Purpose 

The  purpose  of  this  equipment  upgrade  was  to  allow  the  TekScan™  sensor  system  to  be 
used  as  an  input  to  the  dynamic  biomechanical  model.  This  was  accomplished  in  two 
ways:  1)  to  write  new  software  to  allow  direct  access  to  the  raw  data  from  the 
TeKScan™  thus  creating  the  potential  to  extract  specific  contact  force  histories  from 
Tekscan™  and  2)  to  develop  the  ability  to  determine  normal  force  vectors  and  locate 
them  in  space  on  the  body  by  conducting  a  total  body  mapping  of  the  50th  percentile  male 
mannequin. 

6.3  Methodology 

In  the  past  we  have  had  concerns  with  the  accuracy  of  the  Tekscan™  pressure 
measurement  system  and  found  errors  as  high  as  15%  to  20%  in  certain  conditions,  such 
as  bending  the  Mylar  sheet  sensors  around  the  shoulder.  The  system  uses  electro- 
conductive  ink  that  carries  an  electrical  current.  When  a  force  is  applied  to  the  electro- 
conductive  ink,  its  resistance  changes.  The  Tekscan™  software  is  designed  to  record  this 
voltage  change,  identify  the  sensor  location(s),  convert  the  measured  voltage  into  a  force 
value  and,  using  the  known  area  of  the  sensei,  calculate  a  corresponding  pressure  value. 

Any  variation  in  the  thickness  or  resistive  properties  of  the  ink  layer,  either  because  of 
manufacturing  tolerances  or  due  to  bending  of  the  sensor  will  effect  the  response  of  the 
sensei.  Tekscan™  software  algorithms  for  calibration  and  “equilibration”  attempt  to 
compensate  for  the  sensei  to  sensei  variation  in  response.  There  is  also  an  assumption  that 
the  applied  force  is  perpendicular  to  each  sensei  which  may  not  be  true.  Any  of  these 
assumptions  inherent  in  the  Tekscan™  software  calculations  can  lead  to  errors.  Tekscan 
pressure  outputs  continue  to  be  useful  to  observe  ‘hot  spots’,  evaluate  designs  and 
compare  pressure  measurements  to  human  tolerance  guidelines. 

One  goal  of  the  new  approach  in  DBM  is  to  use  the  Tekscan™  pressure  measurement 
system  outside  the  limitations  of  their  commercial  clinical  software.  Tekscan  has  proven 
to  be  very  useful  in  capturing  contact  pressure  distributions  and  in  detennining  the 
relative  magnitudes  of  these  pressures.  The  shortcoming  of  the  technology  has  been  its 
inability  to  determine  the  absolute  magnitudes  of  the  relative  pressures.  Tekscan  also 
provides  a  high  spatial  resolution  at  low  cost,  permitting  the  user  to  pinpoint  the  location 
of  high  pressures.  If  the  pressure  distribution  profiles  measured  with  Tekscan  are 
augmented  with  the  accurate  force  measurement  of  the  exitsing  AMTII  load  cell  in  the 
load  distribution  mannequin,  interface  forces  can  be  modeled  reliably. 
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6.3.1  Software 


The  Tekscan™  software  is  very  sophisticated  in  handling  the  massive  amounts  of  data 
that  emerge  from  sensors  which  consist  of  96  senseis  that  can  be  collected  at  a  rate  of  200 
Hz.  One  way  Tekscan™  manages  this  process  is  to  provide  the  user  with  mainly  graphic 
output  (although  ASCII  output  is  also  available).  Our  goal  was  to  access  the  data, 
including  the  complete  time  history,  calibrate  it  and  convert  it  into  force  values  so  that  the 
modeling  program  can  use  the  output  from  these  three  dimensional  data  arrays. 

The  best  method  at  this  time  appears  to  the  use  a  combination  of  MatLab  (a  numerical 
computation  and  data  visualization  program)  and  Microsoft  Excel  macros.  This  software- 
based  approach  was  developed  with  strong  support  from  research  assistants  who 
developed  a  series  of  user-friendly  data  transfer  files  that  began  with  receiving  data  from 
the  Tekscan™  system  and  using  Matlab  to  manipulate  the  files  into  the  appropriate 
fonnats  compatible  with  MSC  Visual  Nastran™  4D  and  2D  Working  Model.  To  allow 
checking  of  these  data  and  proper  fonnatting,  Excel  macros  were  also  used.  This  software 
is  the  pre-modeling  engine  that  will  assist  us  in  working  outside  the  Tekscan™  software 
to  create  the  DBM.  The  User’s  Technical  Manual,  entitled  “A  User’s  Guide  to  LCSim 
Pressure  Data  Manipulation”  is  included  as  Appendix  A. 

Excel  macros  are  being  used  to  create  a  table  of  force  values  and  the  location  of  the  sensei 
that  measured  that  force.  This  force  array  will  be  needed  in  two  places:  first,  to  map  the 
location  of  each  senseis  onto  the  LC  Simulator  mannequin  and  seconding  to  detennine  the 
normal  force  vector  that  related  to  each  sensei.  Using  the  location  and  magnitude  of  the 
force  vectors,  the  Tekscan™  system  will  be  used  to  examine  the  interaction  forces 
between  the  pack  and  the  person  in  dynamic  biomechanical  modeling.  The  validation  of 
this  approach  is  the  deliverable  in  Phase  II,  Part  6b3. 

6.3.2  Mapping  of  the  Mannequin 

The  full  body  mapping  of  the  50th  percentile  mannequin  is  necessary  for  several  reasons: 
it  will  allow  us  to  locate  the  senseis  in  3D  space,  and  it  will  allow  the  calculation  of  the 
normal  reaction  force.  However,  before  acquiring  the  coordinates  of  the  50th  percentile 
male  mannequin,  full  knowledge  of  the  type  of  software  that  would  be  needed  for  DBM 
was  required  so  that  the  surface  model  infonnation  is  stored  in  an  appropriate  fonnat. 

The  decision  was  made  to  use  MSC  Software  Virtual  Nastran  4D  modelling  and  the 
compatible  “Working  Model”  for  rapid  2D  analysis  (if  needed).  After  that  decision,  the 
mannequin  was  taken  to  Toronto’s  Applied  Precision  Inc.  and  scanned  so  that  the  full 
body  mapping  files  are  now  accessible  by  Visual  Nastran  4D  for  validation  of  the  waist 
belt  model.  The  resolution  of  this  scan  was  2mm  and  it  is  shown  in  the  Figure  6-1  below. 


Equipment  Upgrades  to  Accommodate  Dynamic 
Biomechanical  Modeling 


35 


Figure  6-1.  LC  Simulator  50th  percentile  male  mannequin  3D  scanned  image. 

6.4  Anticipated  Outcomes 

In  future  stages  of  the  DBM,  a  standardized  lower  torso  model  will  be  used  to  validate  the 
approach.  As  this  is  a  symmetrical  form,  the  results  can  be  derived  both  mathematically 
and  from  a  complete  Tekscan™  force  calculation.  The  hip  belt  sensor  should  capture  the 
resultant  force  value.  For  the  shoulders,  the  individual  sensing  elements  can  be  used  to 
find  the  normal  force  vectors  and  sum  them  to  get  the  resultant  shoulder  force  vector. 

This  approach  can  also  be  used  to  validate  the  static  biomechanical  model’s  shoulder  and 
waist  reaction  forces  using  the  normal  force  from  the  Tekscan®  sensors. 

Tekscan™  has  a  rapid  response  rate  and  hence  will  be  very  useful  for  mapping  the  pack- 
person  interface  characteristics.  This  will  be  developed  as  Phase  III. 
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Appendices 

Appendix  A  -  Moment  of  Inertia  Governing  Equations 

The  following  equation  is  used  to  solve  the  moment  of  inertia: 


I  = 


0.2485r 

7 


(m  T " 


Where  : 

/  =  mass  moment  of  inertia  (£g  •  m2) 
m  =  mass  of  plate/object  combination  (kg) 
T  =  period  of  plate/object  combination  (s) 
/» p  -  mass  of  plate  (kg  ) 

Tp  -  period  of  plate  (s) 
r  -  effective  radius  of  plate  (m) 
l  =  length  of  filaments  (m) 
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Appendix  B  -  Centre  of  Gravity  Location  of  a  Measured  Object 


Figure  4-3.  The  top-view  photos  are 
manipulated  with  drawing  software,  such  as 
CorelDraw.  The  intersection  of  the  red  lines 
indicates  the  centre  of  gravity  of  the  object 
in  that  view. 

If  the  proposed  laser  light  assembly  is  used, 
these  steps  are  not  required.  The  object 
can  be  physically  marked  at  the  point 
where  the  centered  laser  light  shines. 


Equipment  Upgrades  for  DBM 


A-2 


Appendix  C  -  Design  Calculations 


Radius  of  Plate 

The  radius  dimension  was  chosen  after  measuring  various  backpacks  and  surveying  the 
curren  t  researchers  of  the  Queen ’s  Ergonomic  Research  Group  backpack  project. 

The  widest  object  the  plate  can  take  is  equal  to  2/3  of  the  height  of  the  triangle  formed  by 
the  cable  attachmen  t  poin  ts. 

The  effective  diameter  of  the  plate  was  chosen  to  be  88.9cm. 

Therefore,  the  effective  radius  is  44.5cm. 


Calculation  of  Ideal  Cable  Length 


r/L  =  constant 

Where  r  is  the  effective  radius  of  the  plate,  and  L  is  the  length  of  the  wire. 

For  the  current  (proven)  pendulum  in  Jackson  Hall  r/L=(24.9cm)/(192.4cm)  =  0.13 
New  pendulum:  L=r/0.13  =  (44.5cm)/0.13  =  342.9cm 
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Appendix  D  -  Calibration 


Calibration  of  Weigh  Scales 

The  calibration  concurred  with  the  manufacturer’s  error  factor  of  +/-llb. 


Weight  (lb)  Scale  A  Scale  B  Scale  C 

5 

5.5 

5 

5 

6 

5 

5 

5.5 

5 

5 

10 

11 

10 

10 

10 

10 

9.5 

10 

10.5 

10 

20 

20 

20 

19.5 

20 

20 

20 

20 

20 

20 

40 

tared 

41 

40 

40 

40 

40 

39.5 

40.5 

40 

40 

5 

5 

5.5 

5 

5 

5 

5 

5 

5 

5 

10 

10.5 

10 

9.5 

10.5 

10 

10 

10.5 

10 

9.5 

20 

20 

20 

19 

20 

19.5 

19 

20 

20 

19.5 

40 

40 

39 

39 

40 

40 

40 

40 

40 

39.5 

120 

tared 

121 

121 

120 

120.5 

120 

119 

120 

120 

119 

10 

10 

9.5 

10 

9.5 

9.5 

10 

10 

10 

9.5 
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Calibration  of  Pendulum 


The  inertia  of  the  empty  pendulum  was  measured  and  compared  with  a  mathematically 
calculated  value.  The  pendulum  was  then  calibrated  with  a  solid  object,  3 1 .4cm  high, 
38.1cm  wide,  19cm  long.  The  mass  was  17.1kg. 


Constants 

Imperia 

1 

Metric 

mass  p= 

22.5 

lb 

10.2 

kg 

total  mass= 

60 

lb 

27.3 

kg 

T  (plate)= 

2.99 

s 

2.99 

s 

1 

(filaments)= 

135 

inches 

3.43 

m 

r  (plate)= 

17.5 

inches 

0.445 

m 

g= 

9.81 

m/s2 

Empty  Plate 

With  Object 

Time  for  10  rotations 

Time  (s)  for  10 
rotations 

1 

29.97 

Trial  1 

Trial  2 

Trial  3 

2 

29.78 

1 

19.92 

19.78 

19.73 

3 

29.79 

2 

19.95 

19.75 

20.19 

4 

29.91 

3 

19.81 

19.96 

19.81 

5 

29.81 

4 

19.93 

19.88 

19.76 

6 

30.03 

5 

19.79 

19.99 

19.73 

7 

29.84 

6 

19.81 

19.73 

19.83 

8 

30.01 

7 

20.03 

20.09 

19.92 

9 

30 

8 

19.92 

19.77 

19.86 

10 

29.98 

9 

19.86 

19.75 

19.79 

11 

29.96 

10 

19.76 

19.93 

19.85 

12 

29.91 

Total 

198.78 

198.63 

198.47 

358.99 

Avg  Period 
(s) 

1.9878 

1.9863 

1.9847 

29.9158 

Period 

2.99158 

Mom  of 
Inertia 

0.239333892 

0.236999 

0.234511 

kgmA 

2 

I  of  plate  alone  (measured) 

Calculated  I 

0.257 

kgmA2 

1= 

1.30821 

kgmA2 

Calculated  I 

1.01 

kgmA2 
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Statistical  Analysis  of  Measured  Periods  During  Calibration 


Anova:  Two-Factor  Without  Replication 


SUMMA 

Count 

Sum 

Average 

Variance 

RY 

Row  1 

3 

59.43 

19.81 

0.0097 

Row  2 

3 

59.89 

19.96333 

0.048533 

Row  3 

3 

59.58 

19.86 

0.0075 

Row  4 

3 

59.57 

19.85667 

0.007633 

Row  5 

3 

59.51 

19.83667 

0.018533 

Row  6 

3 

59.37 

19.79 

0.0028 

Row  7 

3 

60.04 

20.01333 

0.007433 

Row  8 

3 

59.55 

19.85 

0.0057 

Row  9 

3 

59.4 

19.8 

0.0031 

Row  10 

3 

59.54 

19.84667 

0.007233 

Column 

1 

Column 

10 

198.78 

19.878 

0.007307 

10 

198.63 

19.863 

0.015623 

z 

Column 

3 

10 

198.47 

19.847 

0.018112 

ANOVA 


Source  of 
Variation 

55 

df 

MS 

F 

P-value 

F  crit 

Rows 

0.137853 

9 

0.015317 

1.19082 

0.357854 

2.456282 

Columns 

0.004807 

2 

0.002403 

0.186847 

0.831159 

3.554561 

Error 

0.231527 

18 

0.012863 

Total 

0.374187 

29 

Columnl 


Mean 

19.86266667 

Standard  Error 

0.020738843 

Median 

19.84 

Mode 

19.92 

Standard  Deviation 

0.113591322 

Sample  Variance 

0.012902989 

Kurtosis 

1.01809769 

Skewness 

1.051352803 

Range 

0.46 

Minimum 

19.73 
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Maximum 

Sum 

Count 

Confidence 

Level(95.0%) 


20.19 

595.88 

30 

0.04241572 
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Appendix  D,  Annex  #1:  Calibration  Graphs  of  Miniature  Sensors 
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y  x  calc  y  calc  Diff 
x 


15  0  14.4  -0.3  -0.6 
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Actual  Position  vs.  Measured  Position 


Actual  Position  vs  Measured  Position 
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Position  vs  Error 
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Frequency  of  Error  in  the  Y  Dimension 


M  Frequency 
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Appendix  F  -  User’s  Guide  for  the  LCSim  Pressure  Data  Manipulator 


Developed  by  Sabrina  DiVinceza  and  Michael  Hadrovic 
Department  of  Mathematics  and  Statistics 
2001 


What  is  the  Function? 


The  LCSim  Pressure  Data  Manipulator  is  a  Matlab  Graphical  User  Interface, 
which  is  designed  to  aid  in  the  analysis  of  Tekscan  Data,  collected  from  the  Load 
Carriage  Simulator.  The  GUI  calls  upon  a  toolbox  of  Matlab  .m- files  to 
accomplish  this. 


Where  Is  All  This  Stuff? 


There  are  several  Matlab  functions  that  the  LCSim  Pressure  Data  Manipulator  uses  to 
accomplish  the  analysis  of  the  fscan  pressure  data.  They  are  located  on  APLCS#5  in  the 
following  directory,  ‘c:\sab&mike\sab_mike\’. 

Their  filenames  are: 

averagemax.m 

averagemin.m 

cofmass.m 

dudes_wk  1  write  .m 

peaklocation.m 

plotpeak.m 

senselaverage.m 

senselmax.m 

setfile.m 

thesismenu.m 

Also  used  are  the  two  .m-files  that  were  provided  by  Tekscan  to  load  the  data  into 
Matlab.  Loadfsx.m  (edited  slightly)  and  exstr.m  are  located  in  the  directory 
‘  c  :\sab&mike\Matlabfsx\’ . 

The  Macros  which  format  the  spreadsheets  created  by  the  manipulator  are  located  in  the 
excel  file  ‘thesis_VBmacros.xls’,  which  is  in  the  directory  ‘c:\sab&mike\’. 
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How  Do  I  Use  Matlab? 


‘Installing’  the  Data  Manipulator 

APLCS#5  has  Matlab  on  it,  and  the  path  settings  have  been  changed  so  that  Matlab 
knows  where  all  the  necessary  .m-filcs  are.  The  Data  Manipulator  can  be  used  on  any 
computer  with  Matlab  installed  (it  was  created  initially  for  Matlab  version  5.2),  as 
long  as  the  .m-filcs  listed  above  are  on  that  computer  and  you  tell  Matlab  where  they 
are. 

To  do  this,  open  Matlab,  choose  the  ‘File’  menu,  and  choose  ‘Set  Path...’  to  open  the 
Path  Browser.  In  the  Path  Browser,  click  on  the  ‘Path’  menu  button,  choose  ‘Add  to 
Path. . .  ’  and  enter  the  directory  where  the  .m-files  are  stored.  Once  you’ve  done  this  for 
all  the  directories,  Matlab  will  be  ready  to  go. 


Starting  the  Manipulator 

To  begin  using  the  Data  Manipulator,  just  type  ‘thesis  menu’  in  the  Matlab  command 
window.  This  brings  up  the  GUI  from  which  you  perform  all  the  operations  available  in 
the  suite. 


Features  of  the  GUI 

Please  refer  to  “Graphical  Users  Interface”  shown  in  Figure  6-2 

The  column  of  buttons  on  the  left  side  of  the  GUI  window  is  used  to  execute  the  different 

functions. 

•  There  is  a  grey  box  in  the  middle  of  the  window.  This  is  the  description  box. 

After  clicking  on  one  of  the  buttons,  this  box  will  display  an  explanation  of  the 
function  you  just  used. 

•  Below  that  grey  box,  there  is  another.  This  is  the  system  message  box.  This  box 
will  display  messages  regarding  what  the  manipulator  is  doing.  It  will  let  you  know 
about  things  that  are  happening,  but  aren’t  popping  up  in  front  of  your  face. 

•  The  button  in  the  bottom  right  comer  is  the  mystery  button.  See  if  you  can  figure 
out  what  it  does. 

•  Above  this,  there  are  a  bunch  of  yellow  boxes.  These  display  information  about 
the  data  file  you  currently  have  loaded. 

•  Above  these  yellow  boxes  are  two  input  boxes  where  the  user  can  define  the 
threshold  pressure  and  the  mannequin  frequency.  These  values  are  used  in  various 
functions. 

•  To  the  right  of  the  ‘Plot  Peak’  button,  there  is  a  drop  list.  This  list  contains 
plotting  options  only  used  with  the  ‘Plot  Peak’  function. 

•  Finally,  the  input  box  at  the  top,  to  the  right  of  the  Load  .fsx  Data  button,  is  the 
file  box.  This  is  where  the  user  enters  the  path  of  the  data  file  they  want  to  work  with. 
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Manipulation  of  Load  Carriage  Data  through 

MATLAB 


filename  enlry  box 


Figure  6-2  Graphical  Users  Interface  (GUI)  for  LCSim  Pressure  Data  Manipulator 


CHOOSING  A  FILE  TO  WORK  WITH 

To  begin  using  the  functions,  you  must  first  choose  a  data  file  to  work  with.  Enter  the 
filename  in  the  filename  box  located  in  the  upper  right  of  the  GUI.  Be  sure  to  include  the 
entire  path,  the  file  extension,  and  to  place  single  quotes  around  the  whole  thing. 

Example:  in  the  input  box,  type  ‘c:\sab&mike\data.fsx’  to  use  a  data  file  in  that  location. 

Once  the  name  is  entered,  click  on  the  ‘Load  .fsx  Data’  button.  This  will  import  the  data 
file  into  Matlab  and  make  it  the  active  file.  The  results  of  the  load,  including  whether  it 
was  a  success  or  failure,  will  appear  in  the  Matlab  command  window.  Now  you’re 
ready  to  poke  at  the  data  with  other  functions.  If  you  want  to  work  with  a  different  file, 
just  repeat  these  steps  with  the  new  filename. 


What  About  Using  the  Different  Functions? 
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Assuming  you  were  successful  in  loading  a  data  fde,  before  you  actually  start  analyzing 
the  data  you  must  enter  the  mannequin  frequency.  This  is  done  using  the  frequency  input 
box,  which  is  located  just  above  the  yellow  info  boxes  on  the  right.  It  requires  manual 
entry  of  the  frequency  at  which  the  mannequin  was  cycled  up  and  down  during  the  trial 
from  which  the  data  was  gathered.  The  default  value  in  the  box  is  1.8.  The  value  entered 
will  be  understood  to  be  in  Hz. 

Also,  be  sure  to  pay  attention  to  the  description  box  and  the  system  message  box  when 
running  any  of  the  functions.  They  will  indicate  when  there  is  important  information  in 
the  Matlab  command  window,  as  well  as  the  names  of  spreadsheet  files  that  have  been 
created  by  the  functions. 


Plot  Peak 

Activated  with  the  button  of  the  same  name,  the  Plot  Peak  function  plots  a  graph  of 
pressure  vs.  time  for  senseis  that  attained  the  threshold  pressure  at  some  point  during  the 
trial. 

To  use  this  function  you  must  first  set  the  threshold  pressure.  This  is  done  in  the  input 
box  marked  ‘peak  value’  on  the  right  side  of  the  GUI  window.  The  default  value  in  the 
field  is  20. 

Next,  using  the  drop-menu  of  plotting  options,  choose  the  type  of  plot  you’d  like  to  see. 
Once  you’ve  done  this,  you  can  run  ‘Plot  Peak’  by  clicking  the  button.  While  the 
function  is  displaying  the  data  it  may  also  need  you  to  tell  it  what  to  do  by  pressing 
certain  keys.  It  will  request  this  in  the  Matlab  command  window. 

If  the  function  is  run  with  the  ‘multi-view  single  plots’  option  chosen,  it  will  create  a 
figure  full  of  little  graphs,  each  one  corresponding  to  a  different  sensei.  Choosing  ‘single 
plot  at  a  time’  will  place  the  data  for  each  sensei  in  its  own  figure.  The  function  will 
pause  after  it  creates  each  one  and  to  get  the  next,  you  must  press  <enter>  (notice  you 
were  directed  by  the  system  message  box  to  check  the  command  window).  If  you  choose 
‘multiple  plots  per  graph’,  you  will  get  figures  with  one  graph  on  them,  but  there  will  be 
as  many  as  three  lines  on  the  graph,  each  representing  a  different  sensei.  Remember  to 
check  the  command  window  to  see  if  there  are  more  figures  to  display. 


Average  Max 

Average  Max  finds  the  maximum  pressure  that  was  recorded  in  one  of  these  high-pressure 
senseis  for  each  cycle  of  the  mannequin’s  motion.  That  is,  if  the  mannequin  was  moved 
up  and  down  at  a  frequency  of  1.8  Hz,  then  the  averagemax  function  will  find  the 
maximum  pressure  achieved  in  each  10/i8  of  a  second  interval  and  then  average  all  of  these 
values  over  the  whole  trial. 
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This  function,  like  Plot  Peak,  focuses  on  senseis  that  registered  above  the  threshold,  so  be 
sure  to  set  the  ‘peak  value’  accordingly.  Once  this  is  done,  click  the  Average  Max  button 
to  execute  the  function. 

This  data  will  be  displayed  in  the  command  window  and  exported  to  a  spreadsheet  file. 
The  file  is  created  in  the  same  directory  where  the  data  file  is  located.  It  is  also  named 
after  the  data  file  it  was  derived  from.  For  example,  if  the  data  was  from  the  file 
‘c:\sab&mike\data.fsx’,  then  the  spreadsheet  file  will  be  named 
‘c:\sab&mike\data_avgmx.wkl  ’ . 


Average  Min 

Average  Min  operates  in  the  same  manner  as  Average  Max,  but  with  a  few  slight 
differences.  It  targets  the  same  high-pressure  senseis,  but  instead  it  averages  the 
MINIMUM  pressures  from  each  mannequin  cycle. 

The  output  is  also  the  same,  except  the  spreadsheet  file  is  would  be  named 
‘c:\sab&mike\data_avgmn.wkl  ’ . 


Sensei  Average 

When  you  click  on  the  Sensei  Average  button,  the  sensei  average  function  calculates 
mean  pressure  in  each  sensei,  averaged  over  the  entire  trial.  This  information  is  displayed 
in  the  command  window  and  in  a  figure.  The  array  containing  the  mean  pressures  is  also 
written  to  a  file.  Sticking  with  our  ‘c:\sab&mike\data.fsx’  example,  the  file  would  be 
found  with  the  name  ‘c:\sab&mike\data_snslavg.wkl  ’. 


Sensei  Max 

This  function  acts  just  like  Sensei  Average,  but  it  calculates  the  maximum  pressure 
attained  in  each  sensei  during  the  trial.  The  output  is  the  same,  but  the  spreadsheet 
filename  will  be  of  the  fonn  ‘data  snslmx.wkl’. 


Centre  of  Mass 

This  is  actually  Mike  Eklund’s  (Centre  of  Mass  or  footmov)  function  with  modifications 
to  make  it  work  here.  This  function  calculates  the  location  of  the  centre  of  pressure  over 
the  entire  sensei  array  for  each  sample  taken  in  the  trial.  It  then  creates  a  figure  showing 
all  the  pressures  on  the  array  for  one  frame,  as  well  as  the  centre  of  pressure  (marked  with 
an  X).  To  display  the  next  frame,  the  user  needs  to  hit  <enter>.  To  stop  showing  the 
frames,  the  user  enters  a  0.  There  are  prompts  for  these  keys  in  the  command  window. 
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And  that  sums  up  all  the  actions  that  can  presently  be  perfonned  by  the  Manipulator.  Any 
of  the  figures  generated  by  the  functions  can  be  saved  by  choosing  ‘Save’  from  the  ‘File’ 
menu  located  at  the  top  of  each  figure. 


Creating  Excel  Spreadsheets 


The  spreadsheets  created  by  the  Matlab  functions  are  in  a  format  where  their  meaning  is 
not  apparent  and  their  appearance  is  certainly  not  presentable.  For  each  of  the  functions 
that  output  a  spreadsheet,  there  is  a  Visual  Basic  Macro  that  converts  it  into  a  more  useful 
spreadsheet.  These  are  Macros  are  contained  in  the  ‘thesis_VBmacros.xls’  file. 

So,  when  viewing  one  of  the  spreadsheets,  it  is  advisable  to  do  so  in  Excel,  with  the 
‘thesis  VBmacros’  file  open.  Then,  run  the  macro  that  corresponds  to  the  type  of 
function  that  created  the  spreadsheet. 


Quitting  The  Program 


To  exit  the  program,  please  select  the  ‘Quit’  button  located  in  the  bottom  right  corner  of 
the  GUI. 


Equipment  Upgrades  for  DBM 


A-20 


UNCLASSIFIED 


DOCUMENT  CONTROL  DATA 

(Security  classification  of  the  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  document  is  classified) 


3.  TITLE  (The  complete  document  title  as  indicated  on  the  title  page.  Its  classification  is  indicated  by  the  appropriate  abbreviation  (S,  C,  R,  or  U)  in  parenthesis  at 
the  end  of  the  title) 

Development  of  a  Dynamic  Biomechanical  Model  for  Load  Carriage:  Phase  I  Part  A: 
Equipment  Upgrades  to  Accommodate  Dynamic  Biomechanical  Modeling  (U) 

4.  AUTHORS  (First  name,  middle  initial  and  last  name.  If  military,  show  rank,  e.g.  Maj.  John  E.  Doe.) 

Joan  M.  Stevenson;  Susan  A.  Reid;  J.  Tim  Bryant;  Lindsay  J.  Hadcock;  Evelyn  L.  Morin 

6a  NO.  OF  PAGES 

(Total  containing  information,  including 
Annexes,  Appendices,  etc.) 

74 

7.  DESCRIPTIVE  NOTES  (The  category  of  the  document,  e.g.  technical  report,  technical  note  or  memorandum.  If  appropriate,  enter  the  type  of  document, 
e.g.  interim,  progress,  summary,  annual  or  final.  Give  the  inclusive  dates  when  a  specific  reporting  period  is  covered.) 

Contract  Report 

8.  SPONSORING  ACTIVITY  (The  names  of  the  department  project  office  or  laboratory  sponsoring  the  research  and  development  -  include  address.) 

Sponsoring: 

Tasking: 


1 1 .  DOCUMENT  AVAILABILITY  (Any  limitations  on  the  dissemination  of  the  document,  other  than  those  imposed  by  security  classification.) 

Unlimited  distribution 

1 2.  DOCUMENT  ANNOUNCEMENT  (Any  limitation  to  the  bibliographic  announcement  of  this  document.  This  will  normally  correspond  to  the  Document 
Availability  (11),  However,  when  further  distribution  (beyond  the  audience  specified  in  (1 1)  is  possible,  a  wider  announcement  audience  may  be  selected.)) 

Unlimited  announcement 


9a.  PROJECT  OR  GRANT  NO.  (If  appropriate,  the  applicable 

research  and  development  project  or  grant  under  which  the  document  was 
written.  Please  specify  whether  project  or  grant.) 

12CM03 

10a.  ORIGINATOR'S  DOCUMENT  NUMBER  (The official 

document  number  by  which  the  document  is  identified  by  the  originating 
activity.  This  number  must  be  unique  to  this  document) 

DRDC  Toronto  CR  2005-1 1 5 


9b.  CONTRACT  NO.  (If  appropriate,  the  applicable  number  under  which 
the  document  was  written.) 

W771 1-0-7632-01 /A 


10b.  OTHER  DOCUMENT  NO(s).  (Any  other  numbers  under  which 
may  be  assigned  this  document  either  by  the  originator  or  by  the 
sponsor.) 


6b.  NO.  OF  REFS 

(Total  cited  in  document.) 

21 


5.  DATE  OF  PUBLICATION 

(Month  and  year  of  publication  of  document.) 


August  2005 


2.  SECURITY  CLASSIFICATION 

(Overall  security  classification  of  the  document 
including  special  warning  terms  if  applicable.) 

UNCLASSIFIED 


1.  ORIGINATOR  (The  name  and  address  of  the  organization  preparing  the  document,  Organizations 
for  whom  the  document  was  prepared,  e.g.  Centre  sponsoring  a  contractor's  document,  or  tasking 
agency,  are  entered  in  section  8.) 

Publishing:  DRDC  Toronto 

Performing:  Ergonomics  Research  Group,  School  of  Physical 

and  Health  Education,  Queen's  Universityy, 

Kingston,  Ontarioo  K7L  3N6 

Monitoring: 

Contracting:  DRDC  Toronto 


UNCLASSIFIED 


UNCLASSIFIED 


DOCUMENT  CONTROL  DATA 

(Security  classification  of  the  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  document  is  classified) 

1 3.  ABSTRACT  (A  brief  and  factual  summary  of  the  document.  It  may  also  appear  elsewhere  in  the  body  of  the  document  itself.  It  is  highly  desirable  that  the  abstract 
of  classified  documents  be  unclassified.  Each  paragraph  of  the  abstract  shall  begin  with  an  indication  of  the  security  classification  of  the  information  in  the  paragraph 
(unless  the  document  itself  is  unclassified)  represented  as  (S),  (C),  (R),  or  (U).  It  is  not  necessary  to  include  here  abstracts  in  both  official  languages  unless  the  text  is 
bilingual.) 

(U)  Part  A  of  Phase  I  of  the  contract  was  to  develop  the  instrumentation  so  that  a  dynamic 
biomechanical  model  could  be  developed.  The  specific  objectives  were:  a)  to  develop 
Fastrak™  software  for  relative  pack-person  motion,  b)  construct  smaller  strap  force 
gauges,  c)  build  and  calibrate  a  moment  of  inertia  platform,  d)  modify  the  Load  Carriage 
compliance  tester  to  automate  two  degrees  of  freedom  and  e)  create  a  full  body  mapping 
of  mannequin.  This  report  describes  the  purposes  and  outputs  available  from  the 
pack-person  motion  (Section  2.0),  describes  the  development,  construction,  calibration 
and  protocol  for  use  of  smaller  strap  sensors  (Section  3.0)  and  the  moment  of  inertia 
platform  (Section  4.0),  and  the  development  and  steps  involved  in  modifying  the  LC 
compliance  tester  (Section  5.0)  and  the  mannequin  mapping  (Section  6.0).  For  the  most 
part,  tasks  were  developmental  and  construction-based  with  no  data  analyses,  other  than 
to  confirm  the  accuracy  and  precision  of  the  instrumentation. 

This  report  deals  only  with  Part  A  consisting  of  five  sub-parts  of  the  contract.  Part  B 
involved  changing  the  technical  manuals  based  on  the  upgrades  stated  within  this  report. 
Part  C  of  this  contract  is  under  separate  cover  and  was  to  develop  a  long  range  plan  and 
budget  for  dynamic  biomechanical  modeling.  Part  D  was  to  assist  with  the  NATO  PIFM 
Specialist  Meeting  entitled  “Soldier  Mobility:  Innovations  in  Load  Carriage  System  Design 
and  Evaluation”  held  on  27-29  June  2000.  Part  D  is  described  by  the  NATO  RTO  MP-56 
Technical  Proceedings  Report  entitled:  “Soldier  Mobility:  Innovations  in  Load  Carriage 
System  Design  and  Evaluation”. 

1 4.  KEYWORDS  DESCRIPTORS  or  IDENTIFIERS  (Technically  meaningful  terms  or  short  phrases  that  characterize  a  document  and  could  be  helpful  in 
cataloguing  the  document.  They  should  be  selected  so  that  no  security  classification  is  required.  Identifiers,  such  as  equipment  model  designation,  trade  name, 
military  project  code  name,  geographic  location  may  also  be  included.  If  possible  keywords  should  be  selected  from  a  published  thesaurus,  e.g.  Thesaurus  of 
Engineering  and  Scientific  Terms  (TEST)  and  that  thesaurus  identified.  If  it  is  not  possible  to  select  indexing  terms  which  are  Unclassified,  the  classification  of  each 
should  be  indicated  as  with  the  title.) 


(U)  Load  Carriage;  Dynamic  Biomechanical  Modeling;  Compliance  Tester;  Pressure 
Measurement  System;  Fastrak™;  LC  Simulator;  Pack-person  motion 


